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The  purpose  of  this  meeting  is  to  discuss  the ^eS*“‘  imputations  for  propul- 

various  levels  of  engineering models  which  represent  the 
sion.  For  each  level  of  comphcarion  there  ere  their  capabilities  as  well 

The  focus  of  this  meeting  will  be  01,e^01"‘  ^[methods)  which  can  be  applied  to  CFD 
to  higher  order  moment  closure  sc  icm  J devdo{,ing  one.poi„t  closure  models, 

““TfSa  DNS,  will  be  also  discussed  to  some  extent. 


ThU  m!etinTwra  consist  of  three  sessions  and  will  last  about  one  and  half  days 

This  meeting  wm  first  tw0  sessions  each 

Each  session  will  have  three  or  five  y a comment  presentation  (10  minutes) 

position  presentation  (40  min m es)  wi  osition  presentations  (30  minutes)  and  one 

and  a discussion.  In  session  III,  t e speakers  and  the  discussions  will  be 

The  viewgraphs  of  the  presentations  will  be  collected  to  be  d.stnbuted  ater. 


3.  ORGANIZING  COMMITTEE 

T -H.  Shih  (Chairman) 

J.  L.  Lumley  (Honorary  Chairman) 
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M.  Goldstein 
L.  A.  Povinelli 
E.  Reshotko 
J.  M.  Barton 
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4.  OUTLINE  OF  THE  WORKSHOP 
August  21,  1991  (Wednesday) 

08:00-08:15  am  Registration 

08:15-08:30  am  Welcome  by  L.  Povinelli 

Session  I:  Turbulence  Modeling  in  CFD  and  Algebraic  Closure  models. 
Chairman:  E.  Reshotko 

08:30-09:10  am  B.E.  Launder,  “The  current  status  of  turbulence  modeling  in  CFD  and  its 

future  prospects.” 

09:10-09:20  am  D.M.  Bushnell,  “Comment  paper.” 

09:20-09:40  am  Discussion 

09.40-10.20  am  D.  Wilcox,  The  present  state  and  the  future  direction  of  eddy  viscosity 

models.” 

10:20-10:30  am  P.  Spalart,  “Comment  paper." 

10:30-10:40  am  T.  Coakley,  “Comment  paper.” 

10:40-11:00  am  Discussion 

11.00-11.40  am  D.  Taulbee,  The  present  state  and  future  direction  of  algebraic  Reynolds 

stress  models.” 

11:40-11:50  am  A.O.  Demuren,  “Comment  paper.” 

11:50-12:00  am  Discussion 

12:00-01:30  pm  Lunch  Break 

Session  II:  Second  Order  Closure  and  PDF  Method. 

Chairman:  J.L.  Lumley 

1:30-2:10  pm  T.-H.  Shih,  “The  present  state  and  the  future  direction  of  second  order 

closure  models  for  incompressible  flows.” 

2:10-2:20  pm  J.R.  Ristorcelli,  Jr.,  “Comment  paper.” 

2:20-2.30  pm  C.G.  Speziale,  “Comment  paper.” 

2:30-2:50  pm  Discussion 

2:50-3:30  pm  T.B.  Gatski,  “The  present  state  and  the  future  direction  of  second  order 

closure  models  for  compressible  flows.” 

3:30-3:40  pm  J.  Viegas,  “Comment  paper.” 

3:40-3:50  pm  G.  Huang,  “Comment  paper.” 

3:50-4:10  pm  Discussion 

4:10-4:50  pm  S.  Pope,  “The  present  state  and  the  future  direction  of  pdf  methods.” 
4:50-5:00  pm  E.E.  O’Brien,  “Comment  paper.” 

5:00-5:10  pm  J.Y.  Chen,  “Comment  paper.” 

5:10-5:20  pm  Discussion 

6:30-9:00  pm  Banquet  (Pierre  Radisson  Inn,  Great  Northern  Blvd.) 
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Session  I 


Turbulence  Modeling  in  CFD 


and  Algebraic  Closure  Models 


August  22,  1991  (Thursday) 

Session  III:  Unconventional  Turbulence  Modeling. 

y- > — - 

f "• " r* 

09:30-10:00  am  S.  Orszag,  unsteady 

10:00-10:30  am  RJL  Mankbad., 

10:30-11:20  am  W.K.  George,  J.H.  Ferziger  “Tie  role  of 


11:20-11:40  am  Discussion 
11:40-12:00  Concluding  Remarks. 
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Aim  and  Scope 


u i *- — ; — 

1 T etatui  of  turbulence  modelling 

—■  - 

, , r _,h,t  fyoes  of  problem  canfshould  be 
To  rive  concrete  examples  of  VfV 

ucUed  wilh  differed  levels  of  closute  mod 

Pinico,„  emphasis  oo  apphcsfloos  b - — 

treatments 
Models  considered: 

. (Isotropic)  Eddy  Viscosity  Models  (EVM) 

• Algebraic  Second-Moment  Closures 

. Differential  Second-Moment  Closures  (DSM) 

No  space  to  consider  numerical  strategies  needed  for  non-EVM 
treatments 


Basls 

L J- 4-—-—,  « ln  indication  ot  the  turbulen 


' 

,in  = 0 

# Main  industrial  interest  is  in  applying  turbulence  models  In  conditions 
where 'these  conditions  are  not  satisfied! 

, Seem  to  perform  worst  in  2D  curved  flows  and  where  body  forces 
act  in  direction  of  primary  velocity  gradient 

. compressibility  effects  on  turbulence  not  adequate, y accounted  for 
withPan  eddy-viscosity  stress-strain  relation 

• rssa-sf  -- s—  s 

available 
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Eddy  viscosity  models  - EH: 
Near-Wall  Strategy 


• EVM's  rarely  give  satisfactory  levels  of  uju:  away  from  waU 

vicinity:  if  Reynolds  stresses  are  important  there,  second-moment 
closure  is  needed 

• More  difficult  to  devise  suitable  ASM's/DSM's  for  near-wall  sublayer 

• Hence  most  current  research  on  EVM's  concerned  with  treatment  of 
this  Mlow-ReynoIds-numberM  region 

• Log  laws  are  generally  inadequate  

• cven  a mixing-length  scheme  is  better 

• One-equation  models  for  sublayer  currently  seem  a good 
compromise,  especially  if  used  with  a "floating"  k (the  length-scale 
gradient) 


LeRC  8/91 


Eddy  viscosity  models  - IV: 
Low-Re  k- i models 


• Devised  by  reference  to  2 -dimensional  flows  parallel  to  plane  walls 

• Fairly  satisfactory  in  predicting  laminarization  and  diffusion -controlled 
transition 

• Retura  results  of  uncertain  accuracy  when  used  in  3D.  separated  or 
impinging  flows  or  on  curved  walls 

• Need  for  about  40  nodes  across  sublayer  means  that  computations  at 
this  level  for  3D  flows  are  only  just  feasible 

• Further  development  work  still  required,  guided  by  DNS  data  banks 
(Rodi,  Mansour) 

• At  present  it  is  often  better  to  use  a one-equation  EVM  across 
sublayer  blended  to  a two-equation  model  in  fully  turbulent  region 
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_ g/9lj  Eddy  viscosity  models  V^Apphcat 


Problem:  Flow  through  square  " 
ducTTotating  in  orthogonal  mode  ^ 

Relevant  to:  Internal  cooling  of 
turbine  blades 

Importance  of  both  Coriolis  and 
buoyancy  forces  /K 

Experimental  dam  oi  Warner  « B°- 

Iacovides,  Launder  (1991) 

. -,h  « r 67  x 200  grid  covering  half  cross 

3D  parabolic  code  with  3 
section  and  20  hydraulic  diameters 

. , in  core  matched  to  one-equation  low-Re 
Standard  V < m0?  _ 0 9 in  both  regions 

treatment  across  sublayer.  ae  - OM  m 


LeRC 


Second-Moment  Closure  (SMC) 


I 

, o*t  of  rate  equations  for 
Modelling  level  based  on  aPP^  influential  second  moments  (e.g. 
Reynolds  stresses  and  any  otnci 

beat  fluxes) 

at  ** ,evel 

A modelling  level  intrinsically  better  able  to  cope  with  comple 
than  EVM's 

Approximations  needed  for 

* pressure -strain  correlation,  v^i] 

* Dissipation,  qj  (and  hence  f) 

* Diffusion,  d*,j 
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I LeRC  8/91 1 status  of  SMC  - I The  Basic  Mode)  10 

• A simple  closure  based  on: 

Rotta's  linear  retum-to-isotropy  concept  for  non-linear  part  of 
Vfij* 

"Isotropization-of-production"  concept  for  linear  ("rapid")  carts 
°f  «j.  y 

* Daly-Harlow  generalized  gradient  diffusion  hypothesis  for  d- 

Local  isotropy  foi-  «jj,  JJ’ 

* € equation  used  in  k-€  model 

has  been  extensively  applied  in  2-D  and  3-D  subsonic  flows 

• Performance  nearly  always  superior  to  k-e  EVM  - often  markedly 

• Scheme  now  becoming  available  in  many  commercial  rnH« 

(FLUENT,  FLOW3D,  PHOENIC5,  etc.)  c°mmerc,al  codes 


1 

LeRC  | 
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Status  of  SMC  - n The  Basic  Model  (cont'd) 

11 

• Empirical  extension  of  model  to  Iow-Re  sublayer  has  been 
extensively  applied  by  Shima  (1989)  and  colleagues  to  laminarizing 

• This  model  apparently  does  not  do  well  in  high  M boundary  layers 
however  (Huang,  personal  communication) 

• "Wall  echo"  part  of  ^ performs  quite  incorrectly  in  impinging  flows 

• Performs  rather  poorly  in  free  flows  (round/plane  jet  "anomaly"- 

strong/weak  shear  flow  "paradox")  J ’ 
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Application  of  Basic  Model_  - n RotaUng  13 

- 8/91  ^ne  ,-hannel  flow^Laonder  et  ai  (19871, _J 

: are  absent,  Coriolis  effects  on  stress 

co mponents  C^n o nly  agency  provoking  as^meuic  flow 
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T 1 Application  of  Basic  Model  - IH  Axisymmetric  14 

RC  8/91  P£npinging  Jet,  Craft  & launder  (1991) 
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LeRC  8/91  Algebraic  Stress  Transport  Hypothesis  (ASTH)  16 

• By — suitably  approximating  transport  (convection  and  diffusion)  of 
Ujuj  in  terms  of  k transport,  the  closure  becomes  one  where: 

algebraic  equations  are  solved  for  u^u:  *s 
a differential  equation  is  solved  for  k ^ 

• This  is  what  we  mean  by  an  ASM  closure  I 

• Technique  is  most  powerful  where  transport  terms  are  small  i <■ 

in  wall  flows  *” 

• Most  widely  used  ASTH's  not  coordinate  frame  invariant 

• Properly  invariant  versions  have  been  proposed  (Ahmadi,  ICASE)  but 
do  not  so  far  seem  to  have  been  extensively  tested  in  crucial  flows 


LeRC  8/91  APPlication  °f  ASM:  Circular  sectioned 
S-bend,  Abou  Haidar  et  al  (1991) 


ASM/I-1UN 


"fame  • i rn.cz  mtm,  i 
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Status  of  ASM's 


When  transport  is  small,  useful  reduction  in  computational  effort 
achieved  while  retaining  virtually  same  results  as  DSM 

Nature  of  ASTH  is  then  unimportant 

System  of  equations  is  stiffer  than  when  DSM  is  used:  convergence 
is  often  more  difficult  b 

ASM’s  are  on  their  way  out;  not  worth  developing  new  software  for 
this  level  of  model 
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"New"  generation  SMC's 


19 


Approximation  of  fawj ^ °thl sr  turbulence,  2-component 

states  of .Jurbujeiice-  e&-  1501  H 

turbulence,  — 

Proper  frame  indifference 

Extensive  use  made  of  stress  anisotropy  invariants: 

Aj  ■ ajjajj  ; A 3 * aikakjajj 


aij  " (ui  uj  " 3fiij  ukuk)/k 

Extensive  use  made  of  results  of  direct  numerical  simulations 
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Status  of  new  generation  models: 
I Free  flows 


20 


Extensive  teteing  in  boteOgeneons  news  (iCASE.  Steniosd.  Conte, I. 
LeRC.  UMIST) 

Moderate  testing  in  2D  free  shear  Hows 

UMIST  model  tested  in  2D  recirculating  and  colliding  flows 

Far  greater  width  of  applicability  demonstrated  than  with  basic  model 

Only  known  test  for  swirling  recirculating  flow,  not  fully  successful 


LeRC 
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Application  of  UMIST  model  to 
free  shear  flows 
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Flow 


Expt 


l”1-  •>«'  °-To95n0 

Round  J*l  °'09 

rUne  w»k*  °-0*! 

Plan*  plun*  ,1  ,n 

Mixing  1 ay«r  0.16-0.20 


B as  1 c 
Model 
0.100 
0.105 
0.078 
0.078 
0.16 


Colliding  Round  Jets 


0.6  - 
'o.s  ■ 


New 

Model 

0.110 

0.098 

0.098 

0.118 

0.176 


Basic  Model 

UMIST  model 

\ ooo  Data  Witze  (1974) 
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LeRC  8/91  Present  status  of  new  generation  DSM's  27 

. Substantial  advances  demonstrated  over  earlier  models  for  free  shear 
flows 

# Significant  unresolved  (or  Incompletely  resolved)  * 

modelling  near-waU  turbulence  that  limit  range  of  applicability  of  ] 

available  models 

. Much  remains  to  be  done  in  high-speed  flows:  pr«e«t  sugg«^ 
for  modelling  extra  terms  and/or  physical  processes  see  g y 

be  "quick  fixes" 

• Many  improvements  foreseen  over  next  3-5  years 

• New-generation  closures  now  being  incorporated  into  general  purpose 
3D  solvers 
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LeRC 

8/91 

[ Prospects 

28 

• EVM's,  ASM's  and  DSM's  will  remain  in  use  though  with  steady 

decline  in  importance  of  EVM’s  and  ASM’s  in  favour  of  DSM's 

• Improved  versions  of  low-Re  two -equation  EVM's  should  lead  to 

more  reliable  predictions  of  separated  flows  than  at  present 

• New-generation  DSM  closures  will  soon  (2—3  years)  replace  basic 
model  even  in  commercial  codes 

• Further  refinement  of  sub-models  in  second  moment  closures  can  be 
expected  throughout  this  decade 

• Increasing  attention  to  interfacing  SMC  with  higher  order  approaches 

such  as  LES 

• Increasing  use  of  two -time -scale  schemes  providing  distinct  time 

scales  for  large  and  (fairly)  small  eddies 


LcRC  8/91 


A Reminder 


• Extensive  collaborative  testing/assessment  of  turbulence  models 
currently  underway  - coordinated/organized  by  Professor  P.  Bradshaw 
(Stanford)  (with  a little  help  from  JLL  and  BEL) 

• Outcome  of  that  exercise  will  offer  a more  complete  and  objective 
view  of  state  of  turbulence  modelling  than  is  currently  available 
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Center  for  Modeling  o 
Workshop  on  Engineer 


The  current  status 


/ Turbulence  and  Transition 
- ing  Turbulence  Modeling  - 199* 


N92 


24  516 

Q6  &&- 3 


Comment  on: 
turbulence  modeling  in 


CFD  and  its  future  prospects 


by 


D.  M.  Bushnell 


NASA  Langley  Research  Center 
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THE  MAJOR  PROBLEMS  WITH  RSE  CLOSURES 
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wall  functions  not  suitable  for  separated/3-D  flows,  must 
compute  to  wall 
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Space  Administration 
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POSSIBILITY  OF  SHOCK  INSTABILITY  IN  REACTING 
FLOWS 
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examples  of  imbedded  instabilities  (in  turbulent  flows 

a.  Gortler  (dynamic) 

b.  Karman  shedding 

c.  Intersection  region  horseshoe  vortex  interactions. 


Historically  limited  to  simple  geometry/low  ReN-,  but 
both  boundaries  are  expanding 


National  Aeronautics  and 
Space  Administration 
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imensional  Analysis 


CONCLUDING  REMARKS 
Currently  there  are  only  2 ways  to  “make  numbers 
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fixes  which  (tenuously)  can  extend  its  range  of 
applicability 
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Wall  region  treatments 
Compressibility 


SIMPLEX  VIEW  OF  STATUS  - TURBULENCE  MODELING 
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COMPLICATING  FEATURES  — TURBULENT  FLOWS 
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Low  Reynolds  Number  (transition,  walls,  edges) 
Discrete  Organized  Vorticity  (instability,  B.  C.) 


Shock  Wave  Interaction(s) 


Center  for  Modeling  of  Turbulence  and  gl 

Workshop  on  Engineering  Turbulence  Modeling 


DISCUSSION 


J.  Bardina  (to  B.E.  Launder) 

I just  have  a question  for  Prof.  Launder.  You  dismissed  the  algebraic 
Reynolds  stress  models  real  fast.  To  me  a naturrd  way  to  go  from  two , equ 
tion  models  to  higher  models  is  to  go  through  algebraic  relations  firsL 
you  think  that  the  poor  predictions  are  only  due  to  the  poor  models  and  not 
related  to  the  numerical  instability  issues. 


B.E.  Launder  (reply) 

It’s  a matter  of  taste  and  depends  on  the  problem  you  are  looking  at  It’s 
true  that  if  one  is  thinking  of  it  in  that  relation,  the  idea  of  using  just  th 
ame  two  equation  k - e or  * - ui  or  more  complicated  stress-strain  relation 
got  some  appeal.  But  nobody  would  suggest  ASM  is  an  improvement  of 
physics  over  k - t.  The  problem  you  encounter  in  stiffer  equations  make  i 
an  unattractive  level  to  fall  to.  It’s  beginning  to  pass  I could  be  that  for 
a particular  discrete  set  of  problems  it  would  make  a lot  of  sense  On 
overall  if  you  can  not  use  an  eddy  viscosity  model,  you  should  just  bite 
bullet  and  use  the  Reynolds  stress  model. 


A.K.  Singhal 

I would  like  to  make  a comment  about  the  use  of  nomenclative  You  can 
notice  that  even  the  names  used  for  the  Reynolds  stress  models  by  the  invited 
speakers  were  different.  It  would  be  very  helpful  for  industry  if  modelers 
could  use  the  same  nomenclature. 
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EVOLUTION  OF  EDDY  VISCOSITY  MODELS 
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Neither  condition  is  satisfied  in  a typical  turbulent  shear  flow 
"SOLVE  NOW,  WORRY  LATER"  (J.  Cole) 


TWO  EQUATION  MODELS 
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Specific  Dissipation  Rate 
Dissipation  Time  Scale 
Dissipation  Length  Scale 


EQUATIONS  OF  MOTION 
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k / oo  and 


RELATIONSHIP  BETWEEN  c AND 
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Makes  k-£  model  very  stiff  in  the  sublayer  and 
dictates  the  need  for  viscous  damping  terms 
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Like  the  k-*e  model,  the  k-co  model  suffers  from  the 
so-called  round-jet/plane-jet  anomaly 


INCOMPRESSIBLE  FREE  SHEAR  FLOWS 


MODEL-PREDICTED  BOUNDARY  LAYER  STRUCTURE 
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Sublayer  y ^ v/u 
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(c)  Variation  of  wake  strength  with  pressure  gradient 
The  unmodified  “cross-diffusion”  term  distorts  predicted  defect-layer  structure 


LOVES  SAMUEL-JOUBERT! 
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Viscous  Sublayer  Structure 
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DISSIPATION 


WALL  FUNCTIONS  AND  VISCOUS  DAMPING 
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EFFECTS  OF  COMPRESSIBILITY 
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Sarkar  and  Zeman  modifications  fix  something  that 
wasn't  broken. ..flat-plate  boundary  layer 

Perturbation  analysis  explains  why 


Applications 


58 


Perturbation  Analysis  of  the  Wall  Layer 
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where 
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ed  = t*  F(Mt)  f 
F(Mt)  = [M,2  - Mto2]  H(Mt-Mw) 
= 3/2,  Mt  = 1/4 


Applications 

>w  Compressibility  Term 
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UNSTEADY  BOUNDARY  LAYERS 


V- 

o 

■ mmmmt 

u 

o 

o 

0 

a 

— 

c 

0 

0 

£ 

> 

-Q 

l 

c 

0 

0 

> 

0 

0 

E 

CO 

0 

0 

i— 

> 

0 

_0 

0 

O 

> 

"D 

0 

■PWH 

"D 

0 

C 

"O 

D 

O 

O 

E 

n 

G 0 
CO 


3 

i 


<1>  co 

h-  § 


O 

O 

o 


0 

E 

o 


CO 

0 


■O 

c 

CO 


> O 

"D  +- 

CO  ^ 
0 ^ 

a 

0 

CO 


"o 

c 

0 


0 

CL 

c 

0 

E 

0 

v. 

0 

> 

0 

"D 

> 

JQ 

CO 


0 CO 

E 2 
■;z  o 

0 C 

CL  > 
x a> 
0 DC 


as  x 

m o 


O H 
O U 


X 

o 


V)  H 

W o 


x o 

CL  -< 
O H 
*J  U 

M X 

> V) 


62 


2vft  * 1 80* 


y/s 


INCOMPRESSIBLE  SEPARATION 
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Experiment 
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z-u  snocK-separated  flows  can  be  routinely  done  on 
$3,000  desktop  computer  - 3500  grid  points  ->  3 to 
hours 
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The  k-«  model  warrents  more  attention  from  the  turbulence  community  - it  does 
many  things  well  with  the  same  5 closure  coefficients! 


Workshop  on  engineering  turbulence  modeling 
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Comments  on  position  presentation 
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i iiiiimi  i hi  u 


Samuel-Joubert  flow,  computed  with  k 
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• The  model  does  better  than  CS,  BL  and  k-e  in  adverse  pressure  gradient  (lower  C/,  higher 
II).  On  Samuel-Joubert  the  Cf  is  fine,  but  H is  low  (how  did  they  get  H = 1.39  with 

Re  w 5000  and  almost  no  PG?).  In  transonic  flow  the  model  seems  to  place  the  shock  slightly 
in  front  of  the  experimental  shock. 


Samuel-Joubert  flow,  computed  with  the  new  model 
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Summary  and  Comments 


Wall  bounded  flows : The  k — w model  is  probably  best 
Adverse  pressure  gradient  flows 

Compressible  flat  plate  flows,  esp.  with  heat  transfer 

Roughness  - blowing  - transition 

Simplicity  - no  damping  functions 

Numerical  stability  - eg  leading  edge  start 

But:  Separation  and  reattachment  still  are  problems 
(all  models) 

Free  Shear  Flows : The  k — e is  probably  best 

k — w solutions  depend  on  free  stream  w 

k — e more  corrections  and  improvements  available 

This  presentation:  Compressible- hypersonic  flows  (NASP) 

General  discussion  of  k — model  with  corrections 

Comparison  of  model  prediction  for  flat  plate  flows 
(k  — e,  k — oj,  q — o>) 

Comparison  of  model  predictions  for  a separated  ramp 
flow 
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Generalizations  of  the  k — oj  Model 


p — = fitS  - \pkD  - ft*pu>k  + [{p  + <rkpr)ktj]  j 

at  <3 

p — = ajifitS  - CipkD)  - ppu>*  + [(<7-z.fi  + <ru/i7>j]  j 
dt  k 

2 

pT  = a* pk/v,  S — (uij  4*  «;,»  — ^ijUh,k)ui,ji  D — Uk 

ftaseline  Model 

ft*  = 9/100,  ft  = 3/40,  a = 5/9,  a*  = 1, 

- \ rt  2 

<Jb  — lj  = tTu,  — j,  Oi  — 3 

Low  Re  Model  (Transition) 

a,  a*  = fns(-fa),  <jl  = 5/18 
rifvmpressible  Dissip.  (Compressible  free  shear  layers) 

/3,  ft*  = fns(Vk/c) 

Compression  Mod.,  (Compressible  separation) 

Ci  = 2.4,  (^pl/dt  = 0,  L = >/E/w) 

Algebraic  length  scale  (Reattachment  heat  transfer) 
l±T  = a*p\/kL , L — min(>/k/u , .225y) 

Vnrticitv  length  scale  (Incompressible  separation) 

— <x*p\fkL,  L = rnin(>/k/u>,  .%y/k/\Cj\) 
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»exP  w w*>,exp 


Hypersonic  Cylinder-Flare 
35°  flare  angle  M = 7.05 


experiments 
J-L  k-e  model 
— • baseline 

- L-mod 

- - L-mod  + C-mod 

- L-mod  + S-mod 


experiments 
J-L  k-e  model 

— baseline 

- L-mod 

- - L-mod  + C-mod 

— L-mod  + S-mod 


surface  pressure 


surface  heat  transfer 


210-13 


dxew  JA  dx0,ooMo/AO 


Hypersonic  Cylinder-Flare 
35°  flare  angle  M = 7.05 


s (cm)  210-14 
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Woo«exP  w woo.exp 


Hypersonic  Cylinder-Flare 
35°  flare  angle  M - 7.05 


o experiments 

q-o)  model 
baseline 

jl°j  **" 
1°. 

L-mod 

L-mod  + C-mod 

H 

L-mod  + S-mod 

Ip-! 

Wf  surface  pressure 

'' 

experiments 
q-w  model 

— baseline 

- L-mod 

- - L-mod  + C-mod 

— L-mod  + S-mod 


• NA  ° \ 

/ °/X?S: 

p ii 

if 

? / surface  heat  transfer 
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COMMENT  ON  TWO-EQUATION  MODELS 
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FLAT  PLATE  BOUNDARY  LAYER:  Re 
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FLAT  PLATE  BOUNDARY  LAYER:  Re<>  = 1410 


250.0  500.0  750.0 


CHANNEL  FLOW:  Rer  = 395 


91 


102-0 


CHANNEL  FLOW:  Rer  = 395 
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100.0  200.0  300.0  400.0 


CHANNEL  FLOW:  Rer  = 395 
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DISCUSSION 


D.M.  Bushnell 

I would  like  to  make  a comment  pertaining  to  the  problem  of  numerical 
resolution  and  numerical  fidelity.  People  are  showing  all  kinds  of  results  with 
various  models  without  keeping  track  of  how  well  they  axe  doing  numeric  y. 
There  should  be,  at  some  stage,  in  turbulence  modeling  community  some 
agreement  of  some  ? calculation  with  some  standard  code  and  then  stick 
with  that  kind  of  quality  and  fidelity  the  whole  way  through.  Without  this 
I am  not  sure  what  I am  looking  at  quite  honestly. 


R.  Mankabadi  (to  D.  Wilcox) 


I noticed  Wilcox  in  his  talk  gave  k - e an  incomplete  grade  for  the  case 
of  unsteady  boundary  layers.  You  may  like  to  know  that  Howell  (1980)  and 
Ramaprian  (1983)  used  k-e  model  to  calculate  unsteady  boundary  layer  and 
found  that  k-e  model  could  not  predict  this  if  the  amplitude  of  oscillation 
or  frequency  is  high. 


D.  Wilcox  (to  T.Chitsomboom) 

I have  a comment  about  the  T.  Chitsomboom’s  slide  we  just  saw  for  Re  = 
1410.  I have  grid  independent  solutions  for  that  so  I’ll  be  very  suspect  of 

these  results. 


T.  Chitsomboom  (reply) 

I didn’t  mention,  and  you  didn’t  either,  that  k - tv  model  is  quite  sensitive 
to  a?  boundary  condition.  Difficulty  we  encounter  is  how  to  specify  u at  the 
wall.  Solutions  we  have  shown  are  about  in  the  middle  range;  we  can  get 
some  better  results  than  this  and  also  some  worse. 
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D.  Wilcox  (to  T.  Chitsomboom) 


You  have  to  be  careful  with  the  u>  boundary  condition  because  if  this  is 
messed  up  then  the  boundary  layer  is  all  whacked  up  and  I am  sure  that 
that  is  the  case  in  this  computation. 

M.S.  Anand  (to  D.  Wilcox) 

I like  to  preface  my  question  by  mentioning  that  we  do  pdf  methods  and  for 
this  purpose  we  need  the  time  scale  information.  For  this  purpose  we  solve 
u>  equation  (either  a mean  u equation  or  a stochastic  u>  equation).  In  the 
limited  calculations  I have  done,  I have  not  noticed  the  sensitivity  to  free 
stream  u>  you  talked  about.  I have  done  calculations  of  single  axial  jets  with 
or  without  co-flow;  non-turbulent  and  very  low  turbulence  co-flow.  If  there 
is  sensitivity  could  you  clarify  what  the  sensitivity  is  due  to  and  shouldn’t 
there  be  sensitivity  to  k — e models  too  ? 


D.  Wilcox  (reply) 

No  sensitivity  in  k — e because  it’s  just  not  there.  You  can  vary  freestream 
dissipation  all  over  the  place  and  get  the  same  answer.  In  k — u>  model  I am 
cheating  a little  in  one  regard.  I am  always  doing  the  similarity  solutions  and 
not  marching.  There  are  only  two  values  of  u>  which  satisfy  this.  I vary  these 
a little  bit  in  calculations.  Neither  of  these  values  gives  good  spreading  rate 
but  these  more  or  less  bracket  it.  You  are  not  seeing  the  sensitivity  because 
you  are  calculating  the  free  stream  u>.  If  you  change  the  initial  value  of  uj  to 
start  your  calculation  you’ll  see  it. 


M.S.  Anand  (to  D.  Wilcox) 
That’s  what  I am  doing. 
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D.  Wilcox  (reply) 

I’ll  look  into  numerics.  I believe  my  similarity  solutions  are  really  good. 


W.K.  George 

I would  like  to  focus  our  attention  away  from  the  nitty  gritty  of  solutions  and 
focus  on  when  should  we  expect  eddy  viscosity  models  to  work.  Tennekes 
and  Lumely  remind  you  that  it  was  included  in  their  book  to  show  why 
eddy  viscosity  has  problems.  The  fact  that  you  have  a local  model,  it  can 
not  handle  separation  or  flow  pass  separation.  Although  Reynolds  stress 
models  are  a straightforward  increase  in  complexity  but  it’s  a quantum  leap 

in  adding  physics. 


B.E.  Launder  (to  D.  Wilcox  and  T.  Gatski) 

I would  like  to  clarify  a thing  on  Dave’s  and  Tom’s  talk.  Last  time  I saw  the 
k-u  model  it  seemed  to  me  that  the  second  equation  had  in  it  a supple- 
mentary source  term.  I see  Dave  nodding  and  Tom  saying  no. 

D.  Wilcox  (to  B.E.  Launder) 


There  ain’t  none. 


B.E.  Launder  (to  D.  Wilcox) 


So  it’s  cleaned  up  in  the  current  marketed  version.  It  is  very  interesting  step. 
As  I indicated  in  my  talk  if  for  separated  flows  one  wants  accurate  results  we 
need  to  use  an  extra  source  term  in  dissipation  equation  which  will  remove 
difficulty  with  adverse  pressure  gradients.  It’s  not  a problem  of  hig  aci 
number.  It  is  intrinsically  a problem  of  separated  flow;  we  get  too  big .of  a 
length  scale  near  the  wall  and  too  high  of  heat  transfer  coefficient.  An 
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if  separation  is  provoked  by  a shock  wave,  the  separated  flow  could  be  a 
zero  Mach  number  or  a seven  Mach  number.  Survey  shows  when  one  shifted 
attention  from  adverse  to  favorable  pressure  gradients  k — u did  better  in  the 
adverse  pressure  gradient,  and  k — e predicted  much  better  transition  and 
re-laminarization.  Tom  Coakley  said  if  you  want  to  predict  transition  using 
k—u'2,  model  you  have  to  put  Re  effects  which  are  absent. 

D.  Wilcox  (to  Spalart) 


Several  points  were  made  that  I should  answer  to.  First,  a complete  model 
refers  to  terminology  used  by  ??  a few  years  ago.  It  simply  means  you  can  use 
it  with  out  knowing  anything  whatsoever  about  the  flow,  like  an  appropriate 
mixing  length.  Thus  a two  equation  model  is  about  as  simple  of  a model  as 
you  can  get. 


P.  Spalart 


Why  then  isn’t  a one  equation  model  complete? 


D.  Wilcox  (to  Spalart) 


Because  you  still  need  to  specify  a length  scale. 


You  also  seemed  alarmed  that  I was  using  Clauser’s  data  to  tune  the  k — u 
model.  The  perturbation  method  solution  yields  a similarity  solution  that 
demands  that  /3t  be  constant.  You  must  compare  with  data  where  0t  is 
constant.  Whether  or  not  this  is  a limited  data  set,  it  is  the  only  way  that 
is  formally  consistent  with  the  perturbation  solution.  When  I go  to  a non- 
similar solution,  the  perturbation  solution  results  are  certainly  bourn  out. 
That  is  why  I use  Clauser’s  data.  It’s  dictated  by  a mathematical  necessity. 
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P.  Spalart 


How  many  of  the  moments  of  the  turbulence  satisfy  this  right  kind  of  scaling? 
Scaling  which  is  Ur  equals  a constant? 


D.  Wilcox 


Well,  I don’t  know  but  this  shouldn’t  invalidate  the  analysis. 

You  were  worried  about  whether  the  velocity  profiles  and  the  shear  stress 
data  came  from  the  same  calculation-absolutely!  You’re  only  talking  abou 
one  data  point  out  near  the  edge. 


C.  Speziale  (to  Spalart) 

Are  you  saying  that  there  is  no  destruction  term  if  there’s  no  wall?  Essen 
tially  the  destruction  term  disappears? 

P.  Spalart  (reply) 

Yes. 


C.  Speziale  (to  Spalart) 

How  then  would  you  do  in  a non-equilibrium  shear  flow?  According  to  this, 
the  only  thing  that  the  eddy  viscosity  term  can  do  is  grow.  I was  thinking 
specifically  of  a situation  where  diffusion  effects  are  small. 

P.  Spalart  (reply) 

It  wouldn’t  work  there. 
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OUTLINE 


• The  need  for  algebraic  stress  models 
— deficiencies  of  the  linear  gradient 
model. 


• Classical  algebraic  Reynolds  stress 
models. 


• Nonlinear  stress-strain  relations. 


• Critique  of  the  models  and  some  new 
developments. 
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LINEAR  GRADIENT  MODEL 


U{Uj 


2 \ dxj  dxi  ) 
ut  — Cukl  — k2  / e 

= 0.09 

Works  reasonably  well  for  near-parallel  shear 
flows,  however,  C ^ is  not  a constant. 


Poor  representation  of  the  normal  stresses 
and,  hence,  does  not  work  well  for  many  mul- 
tidimensional flows. 
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DEFICIENCIES  OF  LINEAR  GRADIENT  MODEL 

FOR  u-iUj 


• Poor  Representation  of  Normal  Stresses 

o Homogeneous  Shear  Flow,  U = U(y) 

— — o % 

u2  = v~  = w-  - -k 
_ dU 

uv  = 

dy 


o Linear  Gradient  Model  : v — w — 2k/ 3 
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STAGNATION  STREAMLINE  TURBULENCE 


dk 

U—  = P - € 

dx 


where  : P 


-<9I7  -r dV 


■ZL- 


d: 


V‘ 


— — (u2  — 


dU 


dy  V 


y dx 


k — e model  : u2  = 2k /3  — ist 


v2  = 2k/ 3 — i/t 


du 

dx 

dv 

dy 


du  dv\  du  ( d" 

P = vt[--  — ) — = 2i,t 


dx  dy  J dx 


dx 
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ALGEBRAIC  REYNOLDS  STRESS  MODELS 


f ( UiUj  ? ' k j € j *5^'  ) — 0 

More  accurately  describes  the  anisotropy  of 
the  turbulence. 

Account  for  the  effects  of  : 

Longitudinal  surface  curvature 

Corner  Geometry 

Swirl 

Buoyancy 

Rotation 

Maintains  the  simplicity  of  a two-equation 
(k  — e)  model  calculation 


Basic  Assumption  = Local  Formulation  : 

Reynolds  stress  depends  only  on  local 
conditions  - turbulence  is  anisotropic 
only  if  maintained  by  velocity  gradient. 

2 

uiuj  = -k6lJ  if  S{j  = 0 


no 


DERIVATION  OF 
ALGEBRAIC  STRESS  MODELS 


Classical  ASM 

• Derived  from  modeled  Reynolds  stress 
equation 

Nonlinear  stress-strain  models 

• Derived  from  : 

o Two-point  closure  theories 

o Continuum  mechanics 

o Expansion  of  classical  ASM 
formulation 


in 


CLASSICAL  ASM 


Algebraic  Stress  Model  formulated  from  Mod- 
eled Reynolds  Stress  Equation  (Rodi,  1972) 


DuiUj 

Dt 


dTjjik 

dxi 


+ Pij  + $ij  + 


1)  Convection  and  Diffusion  Neglected 

2)  Convection  - Diffusion  proportioned  to 
that  of  the  kinetic  energy  equation 


DuiUj  dTiji 

Dt  dxi 


UiUj  ( Dk 
~k~  ( ~Dt 


UiUj 


UiUj 


— £ ) — Pij  + $ij  + €ij 
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Algebraic  Stress  Model  Continued 


Rodi  (1976) 
= 

UjU 

Shear  Flow 


w2  = 


—uv  = 


DuiUj  dTiji  p . 

-dT  = ^ + p'1+$13+  11 

^(P-e)  = Pij+$ij+€ij 

k 


Cr1-”(  utuj  - -k8ij  ) - C2(  Pij  - - P$ij  ) 

2 1 — C2  Pij/ 6 ~ 3 P/tfijj  ^ 

3 u + Ci  l + £(P/c-l) 


^ + id-c2) 

5*  - ^ - C2) 
^ - |d  - *) 

tfdU 

11  e dy 


C. 

/c3 

P/e  + Cj/2  - 1 

€2 

c. 

fc3 

P/e  + Ci/ 2 - 1 

e2 

A:3 

P/e  + Ci/2  — 1 

'e2 

2 C1/2-l  + C2P/t 

3(1  ( P/e  + Ci/2  — 1 )2 
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ALGEBRAIC  STRESS  MODEL  CONTINUED 


s.  — 

— 


- Ci-|-(  UiUj  - hzSij ) - C2(  Pij  - \p6,  ) 


„ tT.  2 _ , / <917,  dUj 

- c3(  Da  - -PStJ  )-c<k(  — + — 


c, 


@2  + 8 
11 


c = 8<%  - 2 

3 11 


„ 30 Ci  - 2 

C4  = 

55 


Shear  Flow 


Pope  (1975) 


w2  = 


2 2 4C  kz  f dU' 

3 3 p/e  + Cx/2-1  Cfl  e2  \dy 


v2  = 


2 2 10C+1  A:3 

3^  + 3 P/e  + Ci/2  — 1 C^e2 


dy 


w2  = 


2 2 6C+1  k3 

3 3 P/e  + Ci/2  — 1 Me2 


#2/ 


uv 


= -C 


tfdu 

e dy 


Ci/2  - 1 - | ( 11C2  - 4C  - 1 ) P/e 

(P/e  + Cx/2-l)2 


C = ( 1 + -C' 

11  V 2 2 
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ALGEBRAIC  STRESS  MODEL 


Comparison  with  data  for  Homogeneous  Shear  Flow 

Data  from  Harris,  Graham,  and  Corrsin  (1977) 

Neglect  Diffusion  : 

DuiuJ  _ uiu]  Dk 
Dt~  ~ ”T~  Dt 


dutu  n 

u— H- 

ax 

UiUj  dk 

k Udx 

u2 

1.96 

2.09 

V2 

0.88 

0.83 

w2 

1.32 

1.24 

uv 

0.62 

0.62 

Stress  Values  : 

Ci  = 1-8 

Cl  = 1.5 

C\  = 1-8 

C2  = 0.6 

o 

II 

<N 

o 

C'2  = 0.54 

Exp. 

u2  / k 

1.195 

1.050 

0.965 

1.004 

v2  j k 

0.402 

0.366 

0.372 

0.398 

w2/k 

0.402 

0.580 

0.662 

0.598 

—uv/k 

0.356 

0.393 

0.305 

0.298 

G, 

0.068 

0.075 

0.058 

0.057 
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CONVECTION-DIFFUSION  ASSUMPTION 


DuiUj  dTijk  D dij  dTiji 
+ — — =k  — + 


Dt 


dxi 


Dt 


dxi 


'LL {XI  j dT[ 
k dxi 


+ 


UiUj 

~k~ 


where  : a2y 


XL  {'LL  j 2 

~V~3ij 


dTjjk  _ dTi 

dxi  k dxi 


D CL{j 

~Dt 


Then  the  ASM  represents  the  asymptotic  so- 
lution of  the  modeled  Reynolds  stress  equa- 
tion. However,  r — k/e  is  the  local  flow  value 
as  determined  from  the  k and  e transport 
equation. 
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SOLUTION  WITH  ASM  FOR 
MULTIDIMENSIONAL  PROBLEMS 


• Solve  set  of  nonlinear  algebraic  equations  as 
part  of  the  solution  - causes  numerical  prob- 
lems. 

• Simplify  the  ASM  formulation 

o P/4  = 1 

o Neglect  certain  terms 

o Neglect  certain  gradient  terms  for  spe- 
cific problems 


• Expand  into  an  explicit  nonlinear  stress-strain 
form  (NLM) 


UtUj 

k 
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EXPLICIT  STRESS-STRAIN  RELATION 
FROM  ALGEBRAIC  STRESS  MODEL 


UiUj 

~T~ 


(P  — e)  — Pij  + $ij  — 


2 

3 UlJ 


$ij  = - 1 kSij  ) - C2(  Pij  - §P6{j  ) - ■ • • 

P UjUjdUi  uiuj  k dU{ 
e e dxi  k e dxj 


Expansion 


UiUj 

= 


2 c , , „(2)  2 , 

3 C-7  + aij  T + aij  T + 


Ahmadi  (1988) 
Taulbee  (1989) 


Rubinstein  & Barten  (1990) 
Horiuti  (1990) 


UiUj 


- 2C£Sij 


a and  j3  are  directly  determined  from  the  constants  of  the 
RSM.  However,  if  only  terms  through  second  order  are  retained, 
the  expanded  formulation  does  not  accurately  represent  the 
original  ASM.  12q 


NONLINEAR  STRESS-STRESS  MODELS 


General  form  (Speziale.  1991) 

&ij  — S{j 

-4  ax  r2(SiiSij-lS2Sij) 

— 4 a2  r2(  Su  Q,i j + Sji  On  ) 

— 4«3  t2(  On  Oij  — ^ 0 Sij  ) 

— 2 a4  t2  D Sij  / Dt 


r = time  scale  (r  = k/e) 

Lumley  (1969) 

Saffman  (1974) 

Yoshizawa  (1984)  DIA  two-point  closure  theory 
Speziale  (1987)  continuum  mechanics 
Ahmadi  (1988)  expansion  of  ASM 
Rubinstein  & Barten  (1990)  RNG  theory 
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NLM  COEFFICIENT  VALUES 


ax 

OL  2 

<*3 

Nisiziwa/Yoshizawa 

0.090 

0.0709 

0.0159 

0.0960 

(Channel  Flow) 
Speziale 

0.090 

-0.0138 

0.0138 

0 

(Channel  Flow) 
Rubinstein  & Barten 

0.0845 

-0.0570 

0.0120 

-0.047 

(RNG  Theory) 
Rubinstein  & Barten 

0.090 

-0.0523 

0.0198 

0 

(ASM  Expansion) 
Comparison  with  RSM 

0.090 

-0.0015 

0.012 

0 

Homogeneous  shear  flow  : 

Cn  = -au/{Tddly) 

al  ~ a 3 — — f 

P 

i — 1 
h—» 

-1- 

^22  )/ 

CM 

a2  = -K«n  ■ 

-«22)/(<)2 

Irrotational  strain  flow  : 


Cf,  = i(a22  - ou  )/ 

al  = ~ |(  all  + a22  )/  (T^)2 
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STRESS-STRAIN  RELATION  FROM 


EXPANSION  OF  REYNOLDS  STRESS  MODEL 


Reynolds  Stress  Equation  : 


DuiUj  _ _ dT^  + p. . + # . . _ |£  ^ . 


Dt 


dxi 


$ij  = - Ci~~  ( UiUj  - \k8ij)  - a ( Pij  - \P8ij) 


- (3  ( Pij  - | P 8ij)  - 2 7 k Sij 


C2+ 8 „ 8<72  — 2 30<72  — 2 

<*  = — — , P = — 77 — > 7 = 


11 


11 


55 


Dissipation  equation  : 


+ = Ycip-o,/ 

Dt  dxi  1 e 2 e 
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EXPANSION  Continued 


(1  a — P)(auS\j  + djiSu  — | akiSikSij) 

(l  cx  “I-  n^lij  H-  djiQn) 


1)  - (C€l 


Expansion  : an  = a- ] V + a^2V2  + • • • 

J IJ  l] 

Requires  r {SkiSu c)1/2  be  small 

Transport  terms  lead  to  small  second  order 
contributions  and  are  neglected. 
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EXPANSION  Continued 


= a-]V  + + 2 af^rf  + • • • 

,P 

where  t = ( CC2  — 1)  — (C€l  — 1)— 


- 2 CpTSij  - 4 ajr2(  SikSkj  - ^S28ij ) 

— 4a2T2(  Sik£lkj  + SjfeOfcj' ) — 2c*4 t Sij 


4/15 

Ci  + C,2  - 2 + (2  - C£l)P/e 

(l-Q-^)C^ 

” Cl  + 2C£3  - 3 + (3  - 2 Cei)P/c 

(1  — o.  + (3)  

2 JC[  + 2C£3  - 3 + (3  - 2 C£l)P/e) 

Cj* 

2(Ci  + 2 C€3  - 3 + (3  - 2 C€l)P/e) 


EXPANSION  Continued 


• Same  general  form  as  nonlinear  stress-strain 
relation. 


• Model  coefficients  are  determined  from 
Reynolds-stress  model  parameters. 


• There  is  no  (ilnSlij  — ^Q26{j)  term.  Possibly 
the  model  for  the  rapid  part  of  the  pressure 
strain  term  in  the  RSM  is  incomplete.  Stress- 
strain  models  derived  from  DIA  or  RNG  the- 
ory contain  this  term. 


• Series  is  valid  for  small  r S.  Solution  with 
this  model,  even  with  more  terms,  does  not 
give  good  results  near  the  asymptotic  state 
in  homogeneous  shear  or  irrotational  strain 
flows. 
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SOLUTION  VALID  FOR 
SMALL  AND  LARGE  kS/e 


D / a ij  \ 

Dr  17/  “ 


TT 


a 


(2-C£i)^-(C1+C„-2)^- 


Ac. . 

15°  lJ 


- (I  - a-  /3)(auSij  + CLjiSu  - lakiSlkSij) 


(l  — Oi  -f*  "f"  djl^ll  i) 


For 


D fa 


ij  \ _ 


Dr  V r 


0 


Correct  form  for  small  r 5 


^ + 
_ i 7 

r J 


• Right  side  closely  represents  the  asymptotic 
solution  to  the  Reynolds  stress  equation, 
Daij/Dt  — 0 


Linear  algebraic  equations  (P/e  retained  implicitly) 
solved  by  the  method  given  by  Pope  (1975). 
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ALGEBRAIC  STRESS  MODEL 


IMPROVED: 


C\  + C€l  — 2 + (2  — Cei) 


___3  Q 

atj  — 


(1  — a — P)  ^a, iSij  + djiSii  — -auSikSij , j r 


+ (1  - ot  + p)  (auQij  + ajiQu)  r 


STANDARD: 


1 + - 

c 


a 


ij 


- a-  0)  (anSij  + cijiSi,  - ^aH, 


Sufi, 


+ (1  — a + P)  (anQij  + ajiQu)  r 
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NEW  EXPLICIT  ALGEBRAIC  STRESS  MODEL 

For  Two  Dimensions 


- 26^(1  - a - f3)gr2S2  (§<^3)  - SjfJ 

— 2C^(1  — ct  + /?)  g t ^ (SuShj  T S jibin') 

4g/15 

Cfl=  1 - §(1  - a - /3)25252  - 2(1  - a + /?)V^2 


Ci  + Ce2  - 2 + (2  - CCl)P/e 


Model  coefficients  are  determined  in  terms  of 
the  Reynolds-stress  Model  parameters. 
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CONVECTIVE  EFFECTS 


WRITE  EQUATIONS  IN  TERMS  OF  rS 
where,  r — k/e  and  S — (SkiSikP 


T 


Ddij 

IdT 


D 


&ij 


Dt  \tS > 


+ 


1 DtS 

s~Wa,j 


1 DtS 
S Dt 


tDS 

S~Dt 


D ( aiJ  \ 
Dt  \tS ' 


C\  CC2  — 2 + (2  — Cc  Q h 


tDS  1 
5 Dt . 


dij 


rS,, 

15  J 

+ (1  — a 4-  f3)  r (at(Qij  + ajiQn) 


4/15 

C\  + C02  — 2 + (2  — Cq  ) j + 
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Stagnation  Streamline  Flow 
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SOLUTION  FOR  THREE  DIMENSIONS 
Pope  (1975),  Spencer  & Rivlin  (1959,  1960) 

For  CLij  ( S\j , 12 jj  ) — a ( S,  Sit ) 

There  are  Ten  Independent  Symmetric  Tensors. 

T(I)  = S 

T(2)  = SS7  — S2S 

T(3)  = S2  -I{S2}/3 

T(4)  = S22  — I{S72}/3 

T(5)  = S2S2  - S2S2 

T(6)  = S22S  — SS22  — 2I{Sf22}/3 

T(7)  = S2SS22  - S72SS2 

T(8)  = SS7S2  - S2S1S 

T(9)  = S22S2  + S2S22  - 2I{S202}/3 

t(10)  = $7s2s72  - n2s2n 
10 

Then  : a2  ?-  = Q(n)  y(.n) 

n = l 

Where  ^ can  be  functions  of  the  invarients 
{S2}.  {S72},  {S3},  {S12S},  {S22S2} 
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3-D  SOLUTION  Continued 

[Ci  + C€2  — 2 + (2  — C€l)—]aij  = —js—Sij 

k 

-f  (1  — a + /3)—(anQij  -f  ajiCLu) 

k 

(l  ^ /^)  (&tZ S[3  -f-  CLjiSn  " 3 ^kl^lk^ij ) 

1 — a+/3  = (1  + 7C')/11 

1 — a — /3  = (5  — QC^/ll  — * small,  neglect  term 
10 

««  = E Q(n)  4n) 

Tl=l 

k2 

= — *2 Cfi—Sij  — a2g2-^(  Sik^lkj  — ^ik^kj  ) 

k3 

— a6g3 —(Clik^kiSij  + Sik^tki^ij  ~ | Sft2Sij  ) 

/c4 

Q?g  ^~(  kl^lm^mj  ^ik^kl^l  ) 

, _ 4 9_ 1 ~ 2 [l1  ~ a + P)9t\  

>t~  15  l + |[(l-a  + /3)g^]2n2+  [(l-a  + /3)<7^]4(^P 

8 (1  — a + /3) {1  — 2 [(1  — a + fi2} 

15  1 + i[(i_a  + /3)9A]2a2+[(i-Q  + /3);7^]4(n2)2 

<7  = [Ci  + 2C£2  - 2 + (2  - C£l)P/s]_1 
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CONCLUDING  REMARKS 


For  simplicity  and  numerical  solution  purposes, 
explicit  forms  are  highly  desirable. 

An  Algebraic  Reynolds  Stress  Model  and  its  cor- 
responding Nonlinear  Stress  Model  can  be  for- 
mulated which  closely  reproduces  the  Reynolds 
Stress  Model  solution  as  long  as  : 

1)  The  mean  velocity  field  is  not  rapidly 
changing. 

2)  There  are  no  boundary  conditions  or  im- 
posed flow  conditions  which  give  rise  to 
strong  non-local  effects. 

Coefficients  in  the  stress-strain  models  are  not  in 
general  constant  but  depend  on  the  strain  field 
and  the  time  parameter  of  the  turbulence. 

A good  explicit  stress-strain  relation  which  rep- 
resents the  anisotropy  of  the  turbulence  should 
replace  the  linear  relation  now  used  in  the  k-e 
model  for  practical  applications. 
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ROTATION/SWIRL 
GRAVITATION/BODY  FORCES 
STREAMWISE  CORNERS 


Workshop  on  Engineering  Turbulence  Modeling, 
NASA  LERC  Cleveland,  OH,  Aug.  1991 


Workshop  on  Engineering  Turbulence  Modeling, 
NASA  LERC  Cleveland,  OH,  Aug,  1991 
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Most  Primitive  form  of  ASM 


Workshop  on  Engineering  Turbulence  Modeling, 
NASA  LERC  Cleveland,  OH,  Aug.  1991 
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DISCUSSION 


A.  Demuren 

I just  want  to  comment  on  the  k-e  model.  It  appears  the  reason  it  performs 
so  poorly  is  the  value  of  epsilon  at  the  wall.  A very  simple  fix  is  to  eliminate 
epsilon  at  the  wall  and  use  a simple  mixing  length.  This  works  very  well, 
and  gives  the  right  behavior  in  adverse  pressure  gradient,  back  facing  step 
and  separated  flow,  etc.  It  is  quite  an  easy  fix  for  the  k - e model  and  yields 

decent  results. 


Ronald  So 

A comment  about  the  compressible  calculation  with  a fc-e  model  What  we 
have  found  is  that  if  you  do  the  analysis  correctly,  you  can  actually  predic 
compressible  flow  very  well  up  to  Mach  10.  What  Dave  has  shown  up  there 
about  Cf  vs.  Mach  number  is  not  quite  correct.  You  can  get  the  predic  ion 
of  the  adiabatic  and  cool  wall  cases  very  well.  We  have  used  the  baseline 
model  and  Sarkar’s  correction. 


D.C.  Wilcox  (to  B.E.  Launder) 


When  you  gave  your  talk  this  morning,  you  said  that  ASM  suffered  “frame 
invariance” . Could  you  comment  on  this. 


B.E.  Launder  (reply) 

It  depends  on  what  hypothesis  you  use  to  relate  the  convective  transport 
of  stress  to  the  convective  transport  of  strain.  Work  attributed  to  Ro  1 
shows  that  you  get  a different  answer  if  your  frame  of  reference  is  at  rest  or 
rotating  at  a constant  angular  velocity.  You  can  devise  a scheme  Ahmadi 
and  Speziale  have  done  so,  that  is  frame  invariant  and  Dale  Taulbee  was 
talking  about  these  things.  At  the  end  of  it  though,  you  aren  t going  to  ge 
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a better  model  out  of  it. 


D.C.  Wilcox  (to  B.E.  Launder) 

Are  you  saying  that  if  I have  flow  over  a curved  wall,  if  you  forget  to  include 
the  Coriolis  and  centrifugal  forces  as  you  go  over  it,  that  this  is  what  messes 
it  up? 


B.E.  Launder  (reply) 


Yes. 


Something  that  Dale  said  towards  the  end  of  his  talk  he  just  slipped  in  there. 
You  guessed  that  if  you  have  important  diffusive  transport,  then  ASM  won’t 
work.  There  are  many  free  flows  where  diffusive  transport  is  very  important. 
I just  don’t  know  of  a good  algebraic  representation  of  it.  My  feeling  is  that 
if  you  haven’t  already  got  the  software  in  place  for  ASM,  then  you  should 
look  beyond  ASM  for  better  answers. 


D.B.  Taulbee  (reply) 


How  about  all  the  people  who  have  k — e programs  sitting  there.  You  can 
easily  upgrade  them  by  changing  the  explicit  stress-strain  relations.  Not 
everyone  has  access  to  RSM’s.  They  just  can’t  buy  them  because  they’re  too 
expensive.  G.  Huang  (comment  to  D.  Taulbee)  ASM’s  are  just  as  complicated 
to  code  as  RSM’s. 


D.  Bushnell 


Brian  Launder  said  it  very  well:  if  you  have  a situation  that  the  physics  is 
such  that  this  ASM  is  fine  then  it  may  work.  Under  the  NASP  contract,  we 
had  a similar  workshop  about  turbulence  modelling  about  two  and  a half 
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WW  we  asked  was,  “do  you  want  the  wrong  answer  very  easily 
a£ht  answer?"  In  the  NASP  project,  inside  the  scramjet  combustors 
the  flows  are  such  that  we  need  to  go  to  RSM’s  to  get  the  proper  physics. 


W.  K.  George 

Since  Dennis  was  free  to  paraphrase 

colleague  Dale  Taulbee.  If,  for  a hell  of  a 

,Q  to  RSM  and  the  physics  is  bad  for  ASM,  then  tilings  w 
£ a i - « model  and  you  shouldn't  be  us,ng  that  either. 

Also  let  me  add  that  there  is  a lot  of  beating  -^s C IUs  problbly  absolutely 
between  0.98  and  0.95  for  the  spreading  rate  of  jets.  It  is  p y 

impossible  to  determine  this  experimentally. 

Let  me  comment  on 

formulation.  If  one  goes  back  a numbers  like  5.1  are  an  average 

jnir^ 

good.  And  the  theory  used  to  interpret  them  is  not  that  good 
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Basic  Equations  and  Model  Terms  for  Incompressible  Flows 
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Equations  for  the  second  moments: 
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ij  and  €{j  must  be  modeled. 


Heat  flux  equation: 
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. Invariant  — independent  of  coordinate  rotation. 
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. Linearity  principle  of  passive  scalar  turbulence. 


Assumptions 

Turbulence  can  be  characterized  by  energy  containing  scales. 
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Mechanical  and  scalar  dissipation  rates: 


Modeling  of  pressure  correlation  terms: 
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Using  the  solutions  of  above  equations,  we  obtain  (for  homo,  flows) 
for  the  rapid  terms: 
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for  buoyancy  terms: 
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For  slow  term,  following  Lumlev  (1978): 
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Note  that  the  dissipation  tensors  £{j  and  ew  have  been  combined 
into  the  slow  terms. 


Models  for  the  rapid  terms: 
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Reynolds  eddy  structure  model: 
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In  addition  to  Reynolds  stress  equations,  we  need  to  solve  six  extra 
modeled  equations  for  Dij. 


0.2  Sij  + Ci(bikSjk  + bjkSik  ^bijbkiSki) 
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reduce  to  Cornell  model. 


Validation  of  the  models: 
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Two-Dimensional  Channel  Flow: 
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Two-Dimensional  Channel  Flow:  Rer  = 180 


Models  for  the  slow  terms: 


185 


Based  on  UjUj  and  eddy  structure  tensors: 
Reynolds  (1990) 


Lumlev  (1978): 


Shih  and  Lumlev  (1985): 
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Validation  of  the  models: 
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Two-Dimensional  Channel  Flow:  Rer  — 180 


Two-Dimensional  Channel  Flow:  Rer  = 180 


CO 

05 


05 

00 

05 


Pi 

0 

-+J 

o 

Pi 

d 

>> 

o 

3 

33 

0 

^3 

u 

a 

m 

"3 

TJ 

o 


* <s> 

CO 

^1 

II 

cO~  ct> 

v — ✓ ■ c-* 

K 


m 

i> 

05 


5-i 

0 

Pi 

a 3 
5-3 


•r~> 

+ 

3 


CO 

rH|  Ol 


II 


CO  'r> 

V ✓ *t4 

(=3 


m 

t- 

05 


5-t 

0 

"d 

Pi 

pi 

d 

5-3 


0 


P* 

5-3 

^d 

Pi 

d 

Pi 

a 3 

S 

0 

tS) 


co 

t— 

05 


>> 

0 


pi 

5-3 

t3 

Pi 

d 

Pi 
a 5 

S 

0 

t>3 


-o  - 

4 I 

f| 

<t>  I 

QC  I 


5h 

0 

"d 

Pi 

P> 

d 

5-3 

rni 

Pi 

d 

43 

d 

5h 

o 


m 

oo 

05 


0 


pi 

5-3 

"d 

3 

* » — < 
rd 
DO 


05 

00 

05 


5h 
0 
"d 
Pi 
pi 
. d 
5-3 

'd 

Pi 

d 

3—^ 

5-i 

o 


193 


Shill  and  Lumley  (1985): 
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Shih  and  Lumlev  (1985): 
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Validation  of  the  models 
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ssure  temperature-gradient  correlatioi 

(1975)  model; Zeman  and  Lumle; 

lei* Craft  et.  al.  (1990)  model 

d George  (1990) 
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Models  for  turbulent  diffusion  terms: 
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Lumley  (1978): 
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Validation  of  the  models: 
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Two-Dimensional  Channel  Flow: 
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Two-Dimensional  Channel  Flow:  RtT  — 180 


agano  and  Tagawa  (1990): 
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(a)  v’u'v' 


(b)  v’u'2 

Fig.  5 Comparison  of  model  results  for  triple 
velocity  correlations  in  wall  turbulence. 


205 


Fig.  3 Shear  stress  profiles  in  2-D  mixing  layer.  A: 
Bradshaw,  et  A:  Gutmark  It  Wygnanski^2^,  — : 

Present  model. 


Fig.  2 Mean  velocity  profile  in  2-D  mixing  layer.  Umax'- 
the  free  stream  velocity,  y$:  the  position  where  U = 
\Umax-  0;  Bradshaw,  et  al^19!,  — : Present  model. 


Fig.  4 Normal  stress  profiles  in  2-D  mixing  layer.  0»  A.  V 
CastrJ20^,  0,A,y:  Gutmark  it  Wygnanskt^21!,  The 
lines  represent  the  present  model. 
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Fig.  10  Mean  velocity  U/U's  profile  in  axisymmetric  jet. 
Us  llie  centerline  inean  velocity.  0:  Abbiss  ct  al''4'. 
V;  Wvgnanski  & Fiedler'"5',  — : Present  model. 
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Fig  11  Shear  stress  uv/Uj  profile  in  axisymmetric  jet. 
O Rodit261,  X:  Wygnanskiii  Fiedler!25',  A:  Abb.sset 
nl,  — : Present  model. 
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Fig.  12  Normal  stress  u2/C/j  profile  in  axisymmetric 
jet  Legend  as  in  Fig.l  1. 
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. 13  Normal  stress  v- fU\  profile  in  axisymmetric 
Legend  as  in  Fig.l  1 
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Fig.  14  Normal  stress  w2/U J profile  in  axisymmetric 
jet.  Legend  as  in  Fig.l  I. 
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Fig.  6 Mean  temperature  in  heated  round  jet  {A,  Lockwood  & 
Moneib  (1980);  o , Becker  el  al.  (1967);  — , present  model]. 


Fig.  7 Square-root  of  temperature  variance  in  heated  round  jet  [A 
Lockwood  & Moneib  (1980);  o,  Becker  ct  al.  (1967);  — , presen 
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One  of  the  weakest  situations  in  second-order  closure  is  Wall- 
bounded  turbulent  flows: 
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'*  N9  2-24 

J.R.  RISTORCELLI:  THESIS  TITLE  ^vy/ 

Sacond  order  turbulence  simulation  of  the  p ,(V 
rotating,  buoyant  recirculating  convection 
in  thP  Czochraigki  crystal  melt 


1)  PROBLEM  STATEMENT 

20M  for  rotating  turbulence 

2)  THE  RAPID  PRESSURE  MODELS  TESTED 

LINEAR  NONLINEAR 

1)  The  IP  model  3)  Launder's  nonlinear 

2)  The  SSG  model  4)  Shih  & Lumley's  model 

5)  The  2DMFI  model 

3)  THE  COMPUTATIONAL  PROBLEM 

4)  SHORTCOMINGS  AND  FUTURE  WORK 
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THE  CZOCHRALSKI  PROBLEM 


2 2 
Res  = cosR  q/v  Res  < 1 05  Rec  = o>cRr/v  R©c  < 1 05 


Gr  = gPATRc/v2  Gr  $ 1 o1 2 Re  = Gr^2 

Pr  = v/a  Pr  - .01  Ar  - H/Rc 

Bi  = £aTmRc/kAT  Bi  £ 4 A = AT/Tm 

Ma  = -a, jATRsPr/pv2  Ma  < 1 05 
1 < Ro  = Re/ 1 Rec  | ^ °° 


Sb  = Res/Gr<102 


Re  < 106 

Ar  - 1 . 

A - .05 


H 
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THF  COMPIITATIONAI  PROBLEM 


1 7 nonlinear  coupled  pde:s 

MEAN  [ U,  V,  W,T  1 


<uu>  r2<uv>  <uw> 

" <0u> 

20M 

r2<vw>  r4<vv>  r2<vw> 

r2<0v> 

<uw>  r2<vw>  <ww>  _ 

_ <0w>  _ 

[ e,  ee  ] <000> 


30M 


<009> 
<00Uj> 
<0UjUj> 
<UjUjUk>  - 


17  IMFQUALIT1ES  (at  least) 

<ujuj>,  <00>  , dij  = <uiuj><00>  - <0ui><0uj> 

<UaUy>  / (<u«yUy>"^2<UaUa  ) 

<0Ua>  I (<00><u<xua  ) 
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t emp*rdtur« 


Figure  7.5:  Temperature,  streamfunction,  vorticity  and  angular  mo- 
mentum; Solution  with  Ristorcelli  and  Lumley’s  rapid  and  buoyancy 
pressure  models.  Gr  = 1 09,  Re  = 3.16 104,  Re,  = 0.0,  Rec  = 
0.0,  Ro  = oo,  Ma  — 103,  Pr  = 0.01,  Bi  = 2.5,  R,/Rc  = 0.5,  ar  = 
1.0,  T = 80. 
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Figure  7.10:  Twice  turbulence  energy,  q2 , and  scalar  variance, 

< 00  >.  second  invariant,  II,  and  (i.j=3,3)  component  of  the 
anisotropy  tensor,  633;  Solution  with  Ristorcelli  and  Lumley  s rapid 
and  buoyancy  pressure  models.  Gr  = 109,  Re  = 3.1610  , Re,  = 
0.0,  Rec  = 0.0,  Ro  = 00,  Ma  = 103,  Pr  = 0.01,  Bi  = 
2.5,  R,/Rc  = 0.5,  ar  = 1.0,  T = 80. 
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Figure  7.7:  Twice  turbulence  energy,  q2,  and  scalar  variance,  < 
60  >.  second  invariant,  II,  and  (i,j=3.3)  component  of  the  anisotropy 
tensor,  &33;  Solution  with  the  SSG  linear  rapid  pressure  model  and 
Ristorcelli  and  Lumley’s  buoyancy  pressure  model.  Gr  — 109,  Re  = 
3.16 104,  Re,  = 0.0,  Rec  = 0.0,  Ro  — oo,  Ma  = 103 , Pr  = 
0.01,  Bi  = 2.5,  RjRc  = 0.5,  ar  = 1.0,  T = 70. 
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FUTURE  WORK  / SHORTCOMINGS 


1)  WALL  FUNCTIONS 

Physics:  boundary  layers  and  turbulence 
Computations:  grid  dependence 

resolution  of  flow  field 
Low  Re  20M  Models:  too  many  nodes 

2)  RAPID  ROTATION  / STABLE  STRATIFICATION 

inadequate  parameterization  of  cascade 

3)  WHAT  IS  THE  TIME  DEPENDENCE  ? 

Joint  Realizability  couples  rapid  models 
What  does  the  averaging  mean  ? 

Long  time  scale  "coherent”  structures 

4)  COMPUTATIONAL  STRATEGIES 

FDAs  reflecting  realizability 
Include  mean  quantities  in  30M  eq's 
FDAs  with  accurate  time  evolution 
FDAs  reflecting  turbulent  diffusion 
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CONTRIBUTIONS  OF  THE  THESIS 


EXPECTATIONS  OF  2DMFI  MODELS 

1)  Flows  with  strong  body  forces 

stable  stratification 
rotation 

magnetic  fields 

2)  Environmental  shallow  water  flows 

industrial  effluents 

mixing  between  bodies  of  water 


3)  Quasi  two-dimensional  geophysical  flows 

large  scale  ocean  mixing 
regional  atmospheric  modeling 
(mesoscale  variability) 


4)  Unsteady  flows 

time  scales  > integral  time  scale 

unsteady  separation 

large  scale  "coherent  structures"  ?? 
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CONTRIBUTIONS  OF  THE  THESIS 


0)  THE  2DMFI  MODEL 

Models  Xpri  & Xpsir=  f(  b,  b2,  b3),  f(  b,  b2,  b3,  <0u>) 
Satisfies  Realizability,  Joint  Realizability 
Satisfies  2DMFI 

Off-realizability  corrections  exact 
Off-geostrophy  corrections  exact 
Free  parameters  available 
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COMMENTS  ON  THE  PRESENT 
STATE  OF  SECOND-ORDER 
CLOSURE  MODELS  FOR 
INCOMPRESSIBLE  FLOWS 


Charles  G.  Speziale 

ICASE,  NASA  Langley  Research  Center 
Hampton,  VA  23665 


Workshop  on  Engineering 
Turbulence  Modeling 

NASA  Lewis  Research  Center 


August  21  - 22,  1991 


Second-order  closure  models  account  for  history  and 
nonlocal  effects  of  the  mean  velocity  gradients  on  the 
Reynolds  stress  tensor. 

Turbulent  flows  involving: 

• Body  forces  or  curvature 

• Reynolds  stress  relaxational  effects 

• Counter-gradient  transport 

are  usually  better  described. 
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REYNOLDS  STRESS  TRANSPORT 

EQUATION 


drjj 

dt 


where 


^ + uV%j 
dxk 


Tij  = u\u) , 


£ = 1/ 


du[  du[ 
dxjdxj 


— - — - 

dxkdxk 


Cijk  = Uill'jUk  + p'ui^jk  + V uj$ik 


'i  .7 


K = \ra 
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ISSUES  IN  SECOND-ORDER 
CLOSURE  MODELING 


Models  for  11^: 


Typically,  it  is  assumed  that 


II  ij  — eAij(h)  + KA4ijke(  b) 


duk 

dxe 


where 


b 


A'kk&ij 

Til 


(anisotropy  tensor) 


These  models  have  deficiencies  in  rotating 
homogeneous  turbulent  flows  (Reynolds  1989 
and  Speziale,  Sarkar  and  Gatski  1990). 


For  the  return  to  isotropy  problem  in  a rotating 
frame  (with  angular  velocity  0),  these  models 
predict  that  the  second  and  third  invariants  of 
b{j  are  independent  of  0 in  contradiction  of  DNS 
and  RDT  (Reynolds  1989). 


Fig.  2 Typical  RDT  solution  for  th«  fotattan  of  initial'y 
anisotropic  homogeneous  turbulence  (by  T.S.  Shih). 
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For  rotating  homogeneous  shear  flow  in  the 
unstable  flow  regime,  these  models  predict  that  the 
growth  rate  A of  the  flow  defined  by 

jy  „\t  _ \t 

K ~ e , £ ^ e 

is  symmetric  about  its  most  energetic  value  (Speziale, 
Sarkar  and  Gatski  1990). 


n/s 
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Consequently,  if  the  most  energetic  state  - where 
A = Amax  ~ is  placed  at  Cl/S  = 0.25,  the  model  will 
exhibit  similarity  with  respect  to  the  Richardson 
number  Ri  = -2Cl/S(l  - 2 Cl/S).  This  is  in 
violation  of  RDT  and  LES  results  (Speziale,  Sarkai 
and  Gatski  1990  and  Speziale  and  Mac  Giolla 
Mhuiris  1989). 


LES  (Bardina  et  al  1983) 
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Speziale,  Sarkar  and  Gatski  1990  recently  showed 
that  the  general  model  is  topologically  equivalent  to 
the  quadratic  model 

n ij  - —C\eb{j  + C2e  (i bikbkj  - 

+C$KSij  + C±K(bikSjk  + bjkSik 
2 ^mn ^mn j ) d*  C§K(b{k\Vjk  H-  bjfcWik) 

in  plane  homogeneous  turbulent  flows  where 
15 ^ | ( dui/dxj  + duj/dxi)  and  W ij  — - (duxjdxj 

—duj/dxi).  Based  on  these  ideas,  the  SSG  model 
was  developed. 
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The  SSG  model  yields  only  modest  improvements 
on  the  Launder,  Reece  and  Rodi  model.  Substantial 
improvements  will  only  come  if  I Lj  is  taken  to  be  a 
nonlinear  function  of  the  mean  velocity  gradients. 

Two  possible  approaches  are: 

(1)  The  eddy  structure  model  of  Reynolds  (1990) 

(2)  Tensor  dissipation  models  (Speziale,  Raj  and 
Gatski  1990). 
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NEAR  WALL  MODELS 


We  currently  do  not  know  how  to  properly  inte- 
grate second-order  closure  models  to  a solid  bound- 
ary! The  major  problem  lies  in  the  pressure  strain 
correlation  11^  . The  commonly  used  near  wall  mod- 
els for  II ij  have  two  major  deficiencies: 

(1)  The  ad  hoc  dependence  of  nzj  on  the  unit  normal 
rii  to  the  wall.  This  does  not  allow  for  the  proper 
treatment  of  wall  bounded  flows  with  corners. 


(2)  Asymptotic  consistency  is  satisfied  through  sin- 
gular differential  equations;  for  example 

»£?  = + <*»■> 
for  the  near- wall  behavior  of  T12.  This  can  cause 
problems  in  numerically  recovering  an  asymptoti- 
cally consistent  solution. 
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Entirely  new  approaches  are  needed  for  the 
near  wall  modeling  of  II ij\ 

We  are  at  the  end  of  the  road  for  models  of  the  form 


II  ij  = eAijih)  + KMijkt{  b) 


duk 

dxt 


with  ad  hoc  near  wall  fixes. 
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NEEDED  IMPROVEMENTS 

• Models  for  Uij  that  are  nonlinear  in  the  mean 
velocity  gradients. 

• Entirely  new  methods  for  the  integration  of 
second-order  closures  to  a solid  boundary. 

• Incorporation  of  directional  information  into 
the  turbulence  length  scale  (possibly  via  an 
integral  tensor  length  scale). 
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DISCUSSION 


S.  Sarkar  (to  T.-H.  Shih) 


I have  a question  to  Dr.  Shih  about  the  slow  term  pressure  strain  correlation 
comparisons  he  showed.  It  seemed  to  me  that  Rotta  and  our  model  gave  the 
same  results.  That  was  little  surprising  because  the  linear  term  coefficien 
were  different  in  the  two  models.  On  top  of  that  our  model  had  a nonlinear 

term. 


T.-H.  Shih  (reply) 


The  nonlinear  term  is  very  small.  Linear  term  coefficient  for  LRR  model  is 
1.5  and  for  your  model  is  1.7. 

S.  Sarkar  (to  T.-H.  Shih) 

We  have  a paper  in  Physics  of  Fluids  in  which  we  compare  the  two  models 
and  they  are  completely  different. 


T.-H.  Shih  (reply) 

Your  nonlinear  term  can  also  have  opposite  sign  to  linear  terms  thus  giving 
results  similar  to  LRR. 


C.G.  Speziale  (to  T.-H.  Shih) 

You  refered  to  SSG  model  as  a linear  model.  There  is  a coefficient  which 
goes  as  square  root  of  second  invariant  and  also  a term  which  contains  a 
production  term  multiplying  the  anisotropic  tensor.  In  precise  mathematical 
terms  it  is  a quasi-linear  model. 
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T.-H.  Shih  (reply) 

If  the  coefficient  is  constant  the  model  is  linear.  In  SSG  model  the  coefficient 
is  a function  of  second  invariant  and  a production  term. 

J.L.  Lumley  (to  C.G.  Speziale) 

You  refer  to  your  models  as  being  equivalent  to  all  the  other  models  but  only 
in  the  equilibrium  situation.  These  flows  are  never  in  equilibrium.  Would 
you  like  to  comment  on  that. 


C.G.  Speziale  (reply) 


Question  is  how  drastic  are  the  departures.  Then  there  is  this  issue  of  cali- 
brating the  coefficients.  My  motivation  for  doing  SSG  model  was  that  most 
of  the  calibration  we  do  is  from  homogeneous  plane  flows  near  equilibrium 
state.  Since  all  the  models  are  collapsing  to  this  degenerate  form,  my  feeling 
was  to  calibrate  the  model  at  this  sate  and  see  the  differences.  It  seemed  to 
be  reasonable. 


J.R.  Ristorcelli  (to  C.G.  Speziale) 


I have  been  judging  these  models  from  the  point  of  view  of  computability. 
SSG  model  doesn’t  compute  very  well.  It  does  better  in  rotating  situations 
then  it  does  in  the  non-rotating  situations.  I imagine  it  would  do  well  in  in 
homogeneous  shear  flow  situation  from  which  it  was  calibrated.  For  me  I 
built  the  principles  of  realizability  in  the  computation  and  I can’t  compute 
the  flow  with  SSG  model. 


C.G.  Speziale  (to  J.R.  Ristorcelli) 


What  happens? 
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J.R.  Ristorcelli  (reply) 

I get  correlation  (coefficients)  which  are  larger  than  unity  or  eigenvalues  of 
matrices  going  to  zero. 

C.G.  Speziale  ( to  J.R.  Ristorcelli) 

But  no  problems  with  k or  e. 

J.R.  Ristorcelli  (reply) 

Well  k is  the  sum  of  these  eigenvalues. 

C.G.  Speziale  (to  J.R.  Ristorcelli) 

SSG  model  satisfies  limited  realizability.  It  does  guarantee  positive  k and  e. 
T.B.  Gatski  (to  J.R.  Ristorcelli) 

I did  various  calculations  with  homogeneous  shear  flows  using  some  nonlinear 
models  e.g.  Shih-Lumley  model.  It  was  very  difficult  to  use  in  homogeneous 
shear  flow  because  it  was  very  stiffening. 

J.R.  Ristorcelli  (to  T.B.  Gatski) 

What  do  you  mean  by  stiffening? 

T.B.  Gatski  (reply) 

All  the  equations  for  these  flows  are  ode’s.  You  are  using  Range  Ivutta 
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method  and  you  need  very  very  small  steps.  The  only  dilemma  with  making 
assessments  of  turbulence  models  where  you  have  pde’s  is  that  you  can  not 
be  sure  unambiguously  that  there  are  no  problems  with  the  algorithm. 

E.  Reshotko  (to  J.R.  Ristorcelli) 


How  do  you  know  that  these  flows  are  turbulent 


J.R.  Ristorcelli  (reply) 


A lot  of  experiments  have  been  done  to  support  this  e.g.  at  AT&T.  Also  the 
Gr  = 1012  and  Re  = 106  for  these  flows. 


E.  Reshotko  (to  J.R.  Ristorcelli) 


Will  your  equations  with  all  the  turbulence  terms  would  give  a laminar  so- 
lution? 


J.R.  Ristorcelli  (reply) 

Turbulence  would  decay  indicating  a return  toward  a laminar  state. 

B.E.  Launder  (to  T.-H.  Shih) 

Did  the  channel  flow  rapid  term  comparisons  you  showed  include  the  inho- 
mogeneous part  of  the  rapid  term  or  the  wall  reflection  effects? 

T.-H.  Shih  (reply) 


No  wall  reflection  or  inhomogeneous  effects  were  included.  From  the  com- 
parisons may  be  we  can  see  how  to  include  the  inhomogeneous  effects. 
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B.E.  Launder  (to  J.R.  Ristorcelli) 

Regarding  your  choice  of  linear  or  nonlinear  rapid  term,  it  seems  that  the 
nature  of  inter-linkage  between  the  stress  and  dissipation  equations  is  crucial 
in  determining  if  you  get  a steady  state  or  a periodic  behavior. 

J.R.  Ristorcelli  (reply) 

Everything  was  same  and  just  the  rapid  term  model  was  changed. 

B.E.  Launder  (to  J.R.  Ristorcelli) 

Since  some  of  the  e equations  you  are  using  are  not  the  ones  advocated  by 
the  model  originators  so  what  you  were  seeing  wasn’t  the  effect  of  just  a 
change  in  the  rapid  term. 
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necessary'  in  order  to  effectively  employ  shock  capturing  techniques  in  the  numeric 
solution  algorithms. 
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The  SSG  Pressure-Strain  Rate  Correlation  Model 
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Comparison  of  the  Speziale,  Sarkar  and  Gatski  (SSG)  model  and  Fu,  Launder  and  Tselep- 
idakis  (FLT)  model  with  large-eddy  simulations  (LES)  of  rotating  shear  flow  for  different 
ratios  of  the  rotation  rate  to  shear  (G/S):  (a)  LES  of  Bardina,  Ferziger  and  Reynolds, 
(b)  FLT  model,  and  (c)  SSG  model. 
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CFLT  Model  „ „„  ' ssc  Model 
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Variation  of  Anisotropies  with  Strain  for  Homogeneous  Shear  with  Longitudinal  Curva- 
ture (symbols  are  the  experimental  results  of  Holloway  and  Tavoularis,  1992):  a)  Craft, 
Fu,  Launder  and  Tselepidakis  (CFLT)  Model,  b)  SSG  Model. 
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where  V u,^u"u"  is  the  production,  and  e5  the  solenoidal  dissipation. 
DNS  is  used  to  verify  the  model  for  pd!  and  to  calibrate  it. 


Pressure-Dilatation  as  a Function  of  Nonclimensional  Time:  a)  Decaying  Isotropic  Turbu 
lence,  b)  Homogeneous  Shear  Flow. 
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and  uj[  is  the  fluctuating  vorticity. 
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Cel  = 1.44,  C&  = 1.90,  C£  = 0.15 
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Variation  of  Scaled  Turbulent  Dissipation  Rate  Ratio  as  a Function  of  Nondimensional 
Time  (a)  and  Turbulent  Mach  Number  (b). 
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Supersonic  Shear  Layer 
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SSG  pressure-strain  with 

Sarkar  et  al.  compressible  dissipation 
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Variation  of  Normalized  Growth  Rate  with  Convective  Mach  Number  for 
Supersonic  Shear  Layer. 
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Variation  of  Normalized  Stress  Distributions  with  Convective  Mach  Number  for  Supersoni 
Shear  Layer:  (a)  Normalized  Streamwise  Stress,  (b)  Normalized  Shear  Stress. 
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with  the  destruction  of  dissiaption  term. 


Normalized  Skin-Friction  Distribution;  a)  as  a function  of  Mach  number  for  adiabatic 
wall,  b)  as  a function  of  wall  temperature. 


Models  for  the  Reynolds  Mass  and  Heat  Flux 
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Numerical  Solution  of  Compressible  Turbulent  Transport  Equati 

• Discretization  Requirements 
— Accuracy 
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♦ The  extension  to  complex  geometries  can  be  accomplished  in  many  configurations  by 
using  a multiblock,  finite  volume  method. 
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effectively  compared. 
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Outline 


• Viegas 

• Opening  comments 

• An  alternative  RSE  model 


• Results  from  VHVRK  (FRAME)  and  LS  models 

• On  s Hear  layers  — spread  rate  comparisons 

• Results  of  compressibility  corrections 

• Huang 

• General  comments  on  RSE  models 

• compressibility  corrections 

• On  the  law  of  the  wall 

• Concluding  remarks 
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An  alternative  RSE  model: 


FRAME  model  (VHVRK,1983) 


Developed  in  fo^co^reo'ible^turbulent 

'bounSry'ldyeia'whic^i^capable  of  predicting  shoe*  boundary 
layer  flows.) 

Builds  upon  pioneering  incompressible  second  order  model  of  W*. 
Compressibility  effects  included  by: 

I Reintroducing^non-zer^divergence  terms  that  were 

* Accounting  fo/non^tro^s^weighted  fluctuating 
vpiocitv  - compressibility  terms. 


Near  wall  af £act=  ^f^^re-raJe'of  strain  terms 

♦ wal^damping  of  quadratic  return  to  isotropy  (slow)  terms 

* Using" the'patre^waged  form  of  the  Hanjallc-taunder 

modifications"^*  ££  t^s^nf fusing  wail 

H.moino  in  the  destruction  part  of  this  equation. 


* Uses  total  energy  equation  (including  k) 

* Uses  "total”  turbulent  dissipation  rate  equation 


Successfully  applied  to: 


* Adiabatic  flat  plate,  M«3,  to  develope  model 

* Supersonic  expansion  at  M-1.76  Dussauge 

* Transonic  shock-boundary  layer  interaction  at  M-1.36  Delery 

* Corner  flow  at  M-3  - Settles 


Lllating  boundary  layer  on  flat  plate  at  M- 
lbatic  and  nonadiabatic  flat  plate  to  M-8  - 


0 - Spalart  -DNS 

Ka rman— S ho e nhe r r 
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SKIN  FRICTION  ON  A NON -ADIABATIC  FLAT  PLATE  AT  MACH  NO. 

KARMAHH— 3CH0BWHBRR  LAW 


ID 

II 


-ioeo* 
[i  1 i • * 

n J* 

lltff 

noc  + X 
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HTUETA 


ALTERNATIVE  CORRELATIONS  OP  SPREAD  RATE 


IMPACT  PRESSURE  THICKNESS 


S'  - GiM.)  (*■  + vWft )) 

i + VW*)  (<LM) 


G{M,)/G(<S) 

Papamoschou  and 
Roshko 


VORTICITY  THICKNESS 


= (Uy  - U3) /(dU/dy)^ 


= Cu(Mc) 


1 - gjj/gi 
1 + £L/Ui 


C„(Afe)/Cw(0) 


Bogdanoff 


velocity  phorle 


TURBULENCE  MODELS  APPLIED 


EDDY  VISCOSITY  (mass  weighted  variables) 

_ k 2 

fif  — Cft  p - 

TURBULENCE  FIELD  EQUATIONS 


CM  = 0.09 


(pk)  t + (p€.jk)  j = -pi+Dk+E), 


(pe)  e + (purfj  = “ C*Pj  +D.  + E, 

£ £t  5 Extra  Compressibility  Terms 
STANDARD  k - e MODEL 


SARKAR  ET  AL.  (SEHK)  k - e MODEL 


Eh  — -aiMt2pe 


Aft  = — 


Ctl  =:  1.0 


Eh=0  and  Et= 0 
Et  = Q 

d — \ffBf 


ZEMAN  k - e MODEL 


Ek  = ~Cd  F(Mt)pe 


E*  = Q 


F(Mt)  = 1.0  - espf—  ((Aft  - 0.  ■)/■•)] 
and  F(Af,)=0,  if  M,  <0.1, 


Cd  = 0.75 


Aft  = 


a = \j2yRT/{l  + 7) 


283 


0/G(0) 


EFFECT  OF  COMPRESSIBILITY  CORRECTIONS  ON 
PREDICTED  FREE- SHEAR  SPREAD  RATES 


0.0  0.5  1.0  1.5  2.0 

Mo 
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EFFECT  OF  COMPRESSIBILITY  CORRECTIONS  ON  IMPACT- 
PRESSURE  THICKNESS  GROWTH  RATE 


o 


(0)0/0 
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Launder  el  ol.  k - e paroboHc  calc 
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Evolution  of  Turbulonco  Modols 


Modeling 


Dpuru'j 

Dt 


dij  + Pij  + $ij  — ptij  + Compressibility  terms 


Incompressible  Modeling 
Turbulence  diffusion,  dj1- 


<£  + 5J5^K)] 


Bn1,1 1 


* g ra.  L _l  gtt h_ t*(-/  i &uju" >1 

$«fc  (.3C«^  e ( 0xm  ^ £»m  flctn  )] 

• Pressure  Interaction,  + $t;,2  + 

• Slow  Term,  #tj,i 


• -cipebij 

• —cipebij  + c^pe^ikbkj  - jiffy) 

• Fast  Term,  $t;,2 

• Launder,  Reece  and  Rodi,  1975 

• Fu,  Launder  and  Tselepidafcis,  1987 

• Craft,  Fu,  Launder  and  Tselepidakis,  1991 

• Sbih  and  Lumley,  1985 


• Speziale,  Sarkar  and  Gatski,  1990 


I 

PE 

i 

K 
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Incompressible  Modeling  (continue) 


• Wall-Reflection  Term,  $ *j> 

• Hanjalic  and  Launder,  1976 

• Gibson  and  Launder,  1978 
♦ Dissipation  rate,  et;- 

• isotropic,  = f with  constants  in  functions  of 
II  and  III  (Launder  and  Shima,  1989  and  Launder, 

1990) 

• anisotropic 

• €i j = [(1  - /)f  Sij  + [Hanjalic  and 

Launder,  1976] 

• models  satisfying  asymptotic  near-wall  behavior 
[Launder  and  Reynolds,  1983;  So,  1991  and  Shih, 

1991] 

• transport  models  model  for  e»j[Kollman,  1991] 

• Heat  Fluxes,  u"0 

• Ui$  ~ ~ 3CT“  0Zh 
—a  . kZTJTTli  &JL 

• UiO  = -CT-UiU'Q^ 

• ASM  type  heat-flux  equation 

• Transport  heat-flux  equations 
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Compressible  Modeling 


• Pressure  dilatation,  p'u jj  ^ 

• Gatski  et  ah,  1991  ( fi(Mt)Pk  and  /2(M()e) 

• Zeman,  1991  (p,2-equation) 

• Rubesin,  1990  (p‘ 2 -equation  + polytropic  process) 

• Taulbee  and  VanOsdol,  1991  (p,2-equation  + modeled 
p'Ufr  Poisson  solution) 

• Fluctuation  velocity  average,  u” 

• Gatski  et.  al.,  1991  (density-gradient  model) 

• Zeman,  1991  (transport  equation  for  u ") 

• Taulbee  and  VanOsdol,  1991  (transport  equation  for 

«'/) 

• Rubesin,  1991  (Constant  total  enthalpy  -f  polytropic 
process) 

• Dilatation  dissipation 

• e = €a  4-  Cd  and  €&  — f(Mt)eM 

• Zeman,  1990 

• Sarkar  et.  al.,  1989 

• Rubesin,  ’’totcil”  e-transport  equation. 

• Rapid  Compression  Model,  pLn  = Constant  (Reynolds, 
1980;  Morel  and  Mansour,  1982;  Voung  and  Coakley,  1987; 
Coakley  and  Huang,  1991;  Rubesin,  1990  and  Zeman, 
1991) 
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Some  Remarks 


• There  are  more  models  than  what  has  been  presented  in 
the  position  paper.  Some  have  been  tested  in  many  "real” 

flows  with  success. 

• Comparison  of  the  models  based  only  on  simple  homogenous- 
type  flows  may  be  misleading. 

• Near-wall  modeling  is  still  a challenging  problem  for  2nd 
moment  closure. 

• Due  to  strong  coupling  among  governing  equations  and 
the  absence  of  numerical  stabilizing  turbulence  viscosity 
in  the  mean-flow  equations,  the  solution  of  Reynolds  stress 
equations  requires  special  attention. 

» Currently,  LRR,  FRAME  and  Launder-Shima  models  have 
been  implemented  in  a N-S  code  and  comparison  of  models 
against  real-flow  experimental  data  is  underway. 

• The  use  the  total  energy  is  necessary  for  hypersonic  flow 
calculations:  Et  ~ E ky  where  E = cvT  4* 


DpE 

Dt 


- Ph  + pe 
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Hypersonic  Cylinder-Flare 
35°  flare  angle  M = 7.05 


surface  pressure 


surface  heat  transfer 


Some  Remarks  (continue) 


• Wall  flows  — The  law  of  the  wall  is  independent  of 
the  Mach  number  if  the  comparison  is  made  based  on 
the  Van  Driest  transformed  variables. 

• Models  using  ^-equation  produce  lower  Von  Kar- 
man  constant } ac. 

• Model  constants  can  be  derived  as  functions  of 
density  gradients, 

• k — u model  is  less  sensitive  to  Mach  number 
effects  (coincidence  ?)  — only  for  wall  flows. 

• Need  more  turbulent  energy  !! 

• Dilatation  dissipation  concepts  make  the  flow  more; 
’’laminar”  — wrong  direction. 

• The  new  pressure  dilatation  model  shown  in  the 
position  paper  also  lowers  ac. 

• Rubesin’s  total  e-model  approach  goes  into  the 
right  direction. 

• Zeman’s  new  pressure  dilation  model  does  an  ex- 
cellent job. 
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Adiabatic  wall,  boundary  layer  flow 


Compressible  law  of  the  wall 
Adiabatic  wall,  couette  flow 
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Some  Remarks  (continue) 


• Modeling  of  the  c equation  is  still  a challenging  problem 
— both  for  incompressible  and  compressible  flows  — Ex- 
perience has  shown  that  a 2-equation-level  model  can  be 
used  to  improve  the  weakness  of  the  e-equation. 

• Two  baseline  test  problems  are  recommended  - one  is  the 
compressible  mixing  layer  and  the  other  is  the  compress- 
ible law  of  the  wall.  Experience  has  shown  that  these  two 
flows  display  completely  different  behavior. 

• Mixing  layer  — As  Mach  number  increases,  the  spread- 
ing rate  decreases. 

• All  unmodified  turbulence  models  fail  to  predict 
this  behavior. 

• This  leads  to  models  designed  to  increase  the  to- 
tal dissipation  rate  as  the  turbulence  Mach  num- 
ber increases. 
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Comment  on  : The  Present  State  and  Future  Direction  of 
Second-Order  Closure  Models  for  Compressible  Flows 


W.  W.  Liou 

Center  for  Modeling  of  Turbulence  and  Transition 
ICOMP/NASA  Lewis  Research  Center 


Turbulent  Dissipation  Rate  - pc 

pc  = p{c3  + cc)  , cs  = vu”u)'{ 

Incompressible  Models  — 0.6  p c3  Uk,k  (in  the  position  paper) 

—r"j  T"j  S'l  + ...  (Liou  and  Shih  (1991)) 

• Preliminary  analysis  (Liou  and  Shih  (1991))  has  shown  that  the  third  order 
moments,  r"  T"-  S",  may  be  as  important  as  the  terms  that  being  retained 
and  may  need  to  be  modeled  as  well. 


Comment  1 : 
D €s  = 

Dt 


Comment  2 : Turbulent  Mass  Flux  - u[ 

Compressibility  Effects:  ..,  reduced  spreading  rate,  ... 

entrainment 


• Turbulent  mass  flux  terms  in  the  Favre-averaged  equations  may  have  a fair 
amount  of  effect  on  the  mean  flow  development  if  they  are  modeled  more  rig- 
orously, especially  for  wall-bounded  flows. 
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Center  for  Modeling  of  Turbulence  and  Transition 
Workshop  on  Engineering  Turbulence  Modeling  - 1991 


DISCUSSION 


B.  E.  Launder  (T.B.  Gatski) 

I must  say  that  I feel  preferably  unpersuaded  by  the  practice  of  adding, 
just  multiplying  the  dissipation  rate  which  comes  out  of  the  dissipation  rate 
transport  equation  by  a factor  that  is  proportions  to  the  Mach  number.  A 
concept  of  the  dissipation  equation  is  that  it  is  really  representing  the  spectral 
transfer  rate  of  energy.  It  is  looking  at  the  large  scales;  the  small  scales  are 
totally  irrelevant.  Admittedly,  you  won’t  get  the  right  behavior  with  the 
so-called  standard  dissipation  equation.  But  surely  one  must  look  at  how  to 
improve  the  transport  equation  rather  than  having  a quick  fix. 

A comment  on  what  John  Viegas  said  earlier.  He  looks  at  two  flows  and 
believes  you  need  a correction  of  different  sign  - wrong!  I believe  we  need 
to  look  at  twenty- two  or  on  hundred  and  twenty- two.  All  we  are  looking  at 
is  the  desperate  sparsity  of  the  compressible  flow  data  base.  That  is  why 
one  can  get  away  with  these  simplistic  ideas.  Just  by  having  et  equal  to  the 
quantity  one  plus  a function  of  Mach  number,  quantity  times  epsilon  won  t 
work. 

S.  Sarkar  (reply) 

First,  this  was  not  meant  to  be  the  only  compressible  fix.  And,  as  far 
as  you’re  saying  that  the  compressible  dissipation  should  not  add  to  the 
solenoidal  dissipation,  because  after  all  our  conventional  wisdom  is  that  this 
is  a fine  scale  thing,  therefore  it  shouldn’t  affect  epsilon.  This  does  not  violate 
that.  What  you  are  saying  is  that  this  has  an  extra  irrotational  component 
which  has  large  scales  and  small  scales.  It  is  not  as  if  the  compressibility  is 
just  changing  the  small  scale  end  of  the  spectrum;  what  it  is  doing  is  creating 
an  irrotational  component  that  is  both  small  and  large  scales.  We  are  just 
choosing  to  look  at  the  small  scales  because  it’s  is  simpler  to  do  it. 


PRECEDING  page  blank  not  filmed 
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Center  for  Modeling  of  Turbulence  and  Transition 
Workshop  on  Engineering  Turbulence  Modeling  - 1991 


G.  Huang  (S.  Sarkar) 

You  can  compose  that  anyway  you  like,  but  only  for  homogeneous  and  con- 
stant property  flow. 

S.  Sarkar  (reply) 

Yes,  absolutely. 
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WORKSHOP  OBJECTIVE 


• "To  discuss  the  present  status  and  the  future 
direction  of  various  levels  of  engineering 
turbulence  modeling  related  to  CFD 
computations  for  propulsion" 

• Combustion  is  an  essential  part  of  propulsion 

• Discuss  PDF  methods  for  turbulent  combustion 


302 


TURBULENT  COMBUSTION  MODELS 


• Essential  to  integrate  the  development  of: 

turbulence  model 
chemical  kinetics 
numerical  method 

• Turbulent/combustion  interactions 

• Tractable  thermochemistry 
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TURBULENT  COMBUSTION  MODELS 
IN  USE  IN  INDUSTRY 


Typically: 


• k-e 

• equilibrium/mixing-limited  combustion 

• finite-volume  codes 
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IMPROVEMENTS  SOUGHT 


• Finite-rate  kinetics  — NOx,  CO,  soot 

— extinction,  ignition 

• Generality  — beyond  idealized  premixed  and 

diffusion  flames 

• PDF  Methods  can  provide  these  improvements 
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PDF  METHODS 


• Solve  modelled  evolution  equation  for  a one- 
point  joint  pdf 

• <j)(x,t)  — compositions,  <|)={<i>i, 

mass  fractions,  enthalpy 

U(x,t)  — velocity 

co(x,t)  — turbulence  frequency  = £/k 

• Hierachy  of  PDF  methods 

- £ 

- U,£ 

U,  co,  ^ 

• Non-linear  reaction  rates  in  closed  form 
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COMPOSITION  JPDF 


• Need  turbulence  model 
(k-£  or  <uiuj>  - e) 

• Gradient-diffusion  model  of  turbulent  transport 
of  (j) 

• e.g.  J.-Y.  Chen  et  al. 

<UiUj>  - £ 

- 4-step  reduced  scheme  for  methane 
solve  for  jpdf  of  5 compositions 
Masri/Bilger/Dibble  piloted  diffusion  flame 


307 


VELOCITY-COMPOSITION  JPDF 


<U>,  <UiUj>  etc.  obtained  from  jpdf 
Need  e equation  (or  equivalent) 


All  convective  transport  in  closed  form  (no 
gradient-diffusion  modelling) 


Connection  to  Reynolds-stress  models 


e.g.  Haworth  & El  Tahry  (GM) 
Anand  et  al.  (Allison  GT) 


FLOW  OVER  A BACKWARD-FACING  STEP 


H - 0.0762  m 1.43 

H3  - 0.254  m Uref  - 0.196  m/s 


MEASUREMENTS:  PRONCHICK  k KLINE  (1983) 

PDF  CALCULATIONS:  ANAND,  POPE  k MONGIA  (1990) 
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y/H 


FLOW  OVER  A BACKWARD-FACING  STEP 
TRIPLE  CORRELATIONS 


-0.71  -0.26  -0.0!  0.47 


-2-1  0 1 -2  -i  0 1 -2-1  0 I -2-1  0 I 

103  <u’2v’>/U^  , 
ret 
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illllllllUl  II 


VELOCITY-FREQUENCY- 
COMPOSITION  JPDF 


• Single,  self-contained  model  equation 

• Describes  distribution  of  £. 

• e.g.  plane  mixing  layer 
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JPDF  CALCULATION  OF  THE  PLANE  MIXING  LAYER 
SCATTER  PLOT:  DISSIPATION  vs.  LATERAL  POSITION 
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PRESENT  STATE  OF  PDF  METHODS 


• Much  research  and  development  work  remains  to 
be  done,  but: 

Realistic  finite-rate  kinetics  have  been 
incorporated 

Applications  have  been  made  to  complex 
2D  and  3D  flows 

Accuracy — should  be  at  least  as  good  as  a 
second-order  closure 
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THE  FUTURE  OF  PDF  METHODS 


• El  Tahry  & Haworth  (General  Motors): 

"...in  our  opinion,  the  PDF  method  is  the  most 
appealing  of  the  one -point  statistical  approaches 
for  in-cylinder  reacting  flows.  Applications  to 
in-cylinder  combustion  can  be  expected  within  a 
few  years. " 


• Correa  (General  Electric): 

"The  prevalent  k-e!  assumed  shape  pdf  closure 
model... must  be  improved  upon  or  replaced 
before  other  quantities  can  be  usefully  predicted. 
An  alternative  is  the  Monte-Carlo/ pdf  approach; 
although  well  proven  for  fully-developed  shear 
flows,  this  method  needs  to  be  adapted  to 
pressure-dominated  flow  in  complex 
geometries. " 


• From  Proposal  for  PDF  research  and 

development  by  Rolls  Royce,  SNECMA, 

MTU to  European  Community: 

...a  joint  velocity-composition  pdf. ..method 
allows  relatively  complex  chemistry  to  be 
simulated  and  also  fully  couples  the  turbulence 
with  the  chemistry.  It  seems  the  only  way 
forward  from  the  present  position. 
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FUTURE  DIRECTIONS  OF 
PDF  METHODS 

1 . Improvements  and  extensions 

2.  Applications  to  practical  combustion  devices 


IMPROVEMENTS 


1 . Reduced  kinetics  — Maas  & Pope  (1991) 


2.  Mapping  closures  — Chen,  Chen  & Kraichnan 

(1989),  Pope  (1991),  Gao 
(1991) 
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EXTENSIONS 


Incorporate  in  jpdf 

Meyers  & O'Brien  1981 

Pope  1991,  Gao  & O'Brien  1991,  Dopazo  1991 


Represent  coupling  between  reaction  and  mixing 


Contains  information  on: 

jpdf  of  £,  and  % (diffusion  flames) 
<c>  and  X (premixed  flames) 


Reconciles  flamelet  and  non-flamelet  approaches 


NUMERICAL  DEVELOPMENT  FOR 
APPLICATIONS 
(WORK  AT  CORNELL) 


1 . Reduced  kinetics — automatic  generation  and 
tabulation  procedures 

2.  Improved  Monte  Carlo/particle  method 

- second-order  accurate  in  space  and  time 

- low  statistical  error 


3 .  General,  robust  pressure  algorithm 


TRANSITION  FROM  k-e  to  PDF 


• Industrial  combustor  codes: 

- complex  geometry — grid  generation 

- models  for  other  processes — sprays,  soot, 
radiation 

- post-processing/integration  in  design 
procedures 


• Incorporate  PDF  methods  within  existing  codes 


• Numerical  method  fundamentally  different 
(particle  method  vs.  finite- volume  method) 
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4 STAGE  TRANSITION 


Starting  point:  finite  volume  code  for 
<U>,  <p>,  k,  8,  <<]>> 

jpdf  of  U,  discard  U information 

(PDF  method  determines  <<])>,  <p>-  incorporate 

reduced  kinetics) 

jpdf  of  U,  <]) 

(finite-volume  code  determines  <p>  and  £) 


jpdf  of  U,  co,  (j) 

(finite-volume  code  determines  <p>) 


jpdf  of  U,  03,  <|> — self-contained  particle  method 


FIRST  STAGE 


finite- volume  Monte  Carlo/ 

code  particle  method 


<U>,  <p>,  k,  8 


<U>,  k,  8 


jpdf  of  (U),  <j> 


• essentially  jpdf  of  (j) 

• but  simple  transition  (2nd  stage)  to 

jpdf  of  <U>,  <|)  (avoids  gradient-diffusion 
modelling) 

• reduced  kinetics  can  be  incorporated 
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CONCLUSIONS 


• Turbulent  combustion  modelling: 
need  to  integrate 

- turbulence  model 

- chemical  kinetics 
numerical  method 


• PDF  methods 

- U,  U,  co,  £ 

- reaction  and  convective  transport  in  closed 
form 

- finite-rate  kinetic  effects 


• Future  model  development: 

mapping  closures 
reduced  kinetics 

add  V<j) 


• Future  numerical  developments: 

more  accurate  particle  methods 
- general  pressure  algorithm 

incorporation  in  combustor  codes 
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COMMENTS  ON 
THE  PRESENT  STATE  AMD 
FUTURE  DIRECTIONS  OE 
F>  D E METHODS 


A J 


E.E.  O'Brien 
SUNY  at  Stony  Brook 


My  first  comment  on  the  presentation  of  S.B.  Pope  is  to  note 
that  Professor  Pope  is  almost  single-handedly  responsible  for  the 
development  of  the  one-point  PDF  method  to  the  state  in  which  it 
can  now  be  reasonably  expected  to  address  actual  engineering 
problems . 

My  second  comment  is  that  I am  in  accord  with  virtually  all 
of  the  points  he  has  made  including  his  first,  which  was  to 
express  surprise  that  a conference  on  "modeling  related  to  CFD 
computations  for  propulsion"  should  be  so  thin  on  combustion 
modeling.  The  PDF  method  he  reviewed  is  relatively  complicated, 
but  it  appears  to  be  the  only  format  available  to  handle  the  non 
linear  stochastic  difficulties  caused  by  typical  reaction 
kinetics.  Turbulence  modeling,  if  it  is  to  play  a central  role 
in  combustion  modeling,  as  it  must,  has  to  be  integrated  with  the 
chemistry  in  a way  which  produces  accurate  numerical  solutions  to 
combustion  problems.  It  is  questionable  whether  the  development 
of  turbulent  models  in  isolation  from  the  peculiar  statistics  of 
reactant  concentrations  is  a fruitful  line  of  development  as  ar 
as  propulsion  is  concerned. 

There  are  three  issues,  two  mentioned  by  S.B.  Pope,  for 
which  I have  prepared  additional  outlines  which  are  appended  to 
this  note. 

a.  The  one-point  PDF  method 

b.  The  amplitude  mapping  closure 

c.  A hybrid  stategy  for  replacing  a full  two-point  PDF 
treatment  of  reacting  flows  by  a single-point  PDF  and  correlation 
(&  cross-correlation)  functions. 

Finally,  I would  like  to  appeal  for  a concerted  effort  to 
obtain  an  adequate  data  base  for  compressible  flow  with  reactions 
for  Mach  numbers  of  unity  or  higher.  DNS  results  have  played  an 
important  role  in  aiding  the  development  of  PDF  models  for 
incompressible  flows.  A similar  role  can  be  played  in  the 
efforts  to  elucidate  the  many  interactions  of  pressure  with  other 
flow  variables  including  species  concentrations. 
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PDF  Method  Outline 


From  1) 

N-S  eqn. 

2) 

Energy  eqn. 

3) 

Equation  of  State 

4) 

Species  conservation  eqn. 

• Generate 

an  evolution  equation  for  the  1 -point  PDF 

(T.S.  Lundgren,  1957;  C.  Dopazo,  1990) 

• Close  the  PDF  equation  where  necessary  by  ’suitable'  closures 


• Use  Monte  Carlo/particle  methods  for  numerical  solution 

(S.  Pope,  1981  +) 

Major  advantage  + linear  increase  of  numerical  effort  with  the 
number  of  dimensions 
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Example:  A 1 -Point  PDFEguation 


single  species  <)>(x,t); 
Statistically  homogeneous  system 


JH 


V<J>  + D V24>+v(<j)) 


1 -point  PDF  equation 

3£(<L_t)_  = _ JL[D£{(V(i))2|$,  p($' 

dt  d$2  ^ 


. DE{  < w : 2 1$}  is  expected  value  of  scalar  dissipation  conditioned  by  the 
scalar  value. 


329 


General  i -Point  PPF_Eguatian 


• Easily  generated,  and  can  include  multispecies  inhomogeneity, 
compressibility,  etc. 

• pfo,$, .b„,p.6;x.  t)  (1) 

• Closed  terms  : advection  and  reaction 


. unclosed  terms  : pressure,  molecular  diffusion  of  all 
quantities  in  (1) 


Closure  strategies  : Satisfy  PDF  realizability  and 

a)  Reproduce  second-order  moment  closures  for  physical  space 
terms  (Pope,  1985) 
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The  Amplitude  Mapping  Closure 


(Chen,  Chen,  Kraichnan  1989;  Gao,  1991) 

• Attractive  for  strongly  non-Gaussian  processes 

• Simplest  example:  Non-reacting  single  scalar  <t>(x,t)  in 

statistically  homogeneous  system. 

• l-Pt  PDF 

3P($,  t)_  = -JL[DE{(V4>)2\$,  t}P($,  t)] 
dC  3$  2 

• Define  6(Z)  time-independent,  homogeneous,  isotropic,  normalized 
Gaussian  r.v.  All  statistics  of  6(z)  are  known  if 

0 = 0,02  = 1 and  fe(r)=0(z)  0(z+r) 


are  given. 

• Define  a scalar  field  4>s(x  ,t)  generated  from  0(z)  from  the  mapping 

x=z/J(t) 

and  <f>s(x,  t)  = X(0(z)  , t) 

• Demand  P($,  t)  = P($s,  t) 


331 


Consequences  of  the  Mapping 


• Since  statistics  of  0(z)  are  completely  known  the  statistics  of  4>s(x,t) 
are  also  completely  known  if  J(t)  and  X are  specified. 


• N & S condition  for  P($,  t)  =P($s,t ) 

Is  £{(V< J>)2|$,  t}  — ( V4>s) 2 |$s,  t) 


• It  turns  out  that  substitution  of  the  mapping  into  the  PDF  equation 
produces  a solvable  equation  for  X 


dx 

3t 


v at>  a<t>2 


J DJ2(t)  ^ 

x is  a normalized  time  scale  dx ri  aL‘ 

Ae 

Note:  J(t)  & *0  (the  on,y  parameter  of  fe  that  matters) 


Appear  only  in  the  time  scale 
/.  the  shape  of  P($,  t)  depends  only 
on  the  mapping  X(0,t). 
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Some  Results  of  the  Mapping 


Symmetric  binary  mixing  (initial  double  delta  PDF) 

P($,0)  = ±6  ($+1)  +±&  ($_!) 
soln:  i?{(V<t>)2|$,T}/F(T)  = exp{-2[er/-1($)  J2} 

Unsymmetric  binary  mixing 

■P(<J>,0  = a6  (<J>+1)  + (l-a)  6 ($-1)  , 0<a<l 

Same  mapping  closure  solution  for  £{(V<f>)  2|$,  t}/f(t) 
General  soln.  has  been  obtained  (Gao,  1991) 


Current  status!  Formal  solutions  have  been  obtained  for 
multispecies  cases  (Gao  & O’Brien,  1991) 

JJyt,  no  reported  success  in  incorporating  it  in  numerical  codes  for 
more  than  one  species 

Also  seems  to  misrepresent  asymptotic  behavior  in  time 


2-Point  PDF 


Advantages 

• Spatial  structures  explicitly  included 

• Self-contained  time  and  length  scales  as  in  spectral 
description  of  turbulence  and  unlike  1 -Point  PDF  & K-e  or 
other  moment  closures 


Disadvantages 

• Dimensions  doubled 

• Closures  harder  to  construct 

• Numerical  work  so  far  limited  to  isothermal  reactions  of  type 
A+B-P. 

Closure  approximations 

Advection  EDQNM 

Diffusion  Linear  Mean  Square  Estimate 

Reaction  None  needed 


New  Wrinkles 

• Use  1 -point  joint  PDF  of  quantities  end  their  gradients 

P(< j>,  V4);x,  t) 


• Hybrid  closures 
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Hybrid  Strategy 


Aimed  at  more  than  2 species,  statistically  homogeneous. 


• Numerical  method  is  fractional  steps  (Yanenko,  1971) 

P(t+At)  - ( I+0AAt)  (I+O^t)  (J+O-At)  P(t) 

I Is  the  Identity  operator 


1 . e. 


■ P°F ' • - 

In  a cycle  of  computations  in  both  composition  & physical  space 

^1?"  has.no  effect  on  1 -Point  PDF  but  it  modifies  the 
elation  and  cross-correlation  functions  f(r,t)  [EDQNM] 

Molecular  diffusion  modifies  the  1 -point  PDF  (LMSE  or 
(knSwnf  SUre’  'f  workab,e  and  the  correlation  functions 

J 2222?  :,eactlon  e,ffects  the  1 -P^nt  PDF  (known)  and, 
similarity)0  m8^  8 corre*a^on  Unctions  (assume 


b) 


c) 


• Reproduces  full  2-point  PDF  results  for 

A + B - P 
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Workshop  on  Engineering  Turbulence  Modeling 


Comments  on  Pdf  Methods 


J.-Y.  Chen 

Combustion  Research  Facility 
Sandia  National  Labs 
Livermore , CA 
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Grand  Challenge  of  Combustion  Engineering 


Challenge:  Significant  Reduction  of  Combustion 

Generated  Pollutants 

Factsi 

■ Pollutant  Formation  Rate  « Fuel 
Oxidation  Rate 

■ Small  Quantity 

■ Highly  Sensitive  to  Interactions 
between  Turbulence  and  Chemistry 

Difficulty: 

■ Not  Capable  of  Solving  Navier-Stokes 
Equations  with  Detailed  Chemistry 

Approach: 

■ "Rational"  Modeling 


i- 
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till  lilill 


Research  of  PDF  Methods  at  Sandia 


• Simple  Geometry  ( Parabolic  Flow) 

• Reduced  Reaction  Mechanisms: 

- Two-step  H2  Flames 

- Three-step  C0/H2  Flames 

- Four-step  CH4  Flames 

- Five-step  CH30H  Flames 

- up  to  six  reactive  scalars 

• Thermal  NO  F ormation  in  Turbulent 
Hydrogen  Jet  Flames 

• Soot  Foramtion  in  C2H4  flames 
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(Fuel  Jet:  45  % CO/15  % H2/40  % N2) 

(Pilot  Jet:  70%CO/30%H2) 

EXHAUST 

RflW 


45  %CO/15%  H2/40  % N2 


Departures  From  Chemical  Equilibrium 


Hydrogen  Methanol 


f 


Carbon  Monoxide  Methane 
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Modeling  Turbulent  Reacting  Flows 


Reynolds  Stress  Model 
for  Velocity  Field 


It  Tt 
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OH  Mo  e Friction  OH  Mo  « F acioi  OH  Mole  Fraction 


Superequilibrium  OH  Radical 


Measurements 


Predictions 


Casa  B,  x/ 0=30 


Case  B.  x/u—30 


Case  A.  x/D-30 


Case  A.  x/D— 30 


: * 

\’t*  ’ V 


Case  A.  x/D—50 


Case  A,  x/D-50 


w.  • 

•Mv- 

vvW  > 

l-V.;  . 

VtC’y*  tL'fc**.’ 

•; f . , . - 


Mixture  Fraction,  f 


Mixture  ^ 


{Mess  Fraction)  YCM4  {Mass  faction!  Y!5U  (Mass  Fraction) 


Scatter  Plot  for  CH4  and  02 


(CH4  Turbulent  Jet  Flames) 


CH4-«ir  FUme  lz 

k/O-20.  c/O— 1 JC3.1.68 
<Wgcf  * nep) 


"k: 


ife  \ 

| « -Y  X. 

iVV'^vfc  -'  - X 

— — ruv  ^ r ’1  ^ i 

o.oo  o.04  o.oe  o.i2  o.ia  0.20  0.24 

(Ma**  Fraction) 
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O 


CH4  air  rUnva  B; 

*/£W0#  #/t)<*f_23,t.4  74l.fl5 
(QUger  4-ifepJ 


5 ®-«- 


CH4  a if  Flame  L: 

*/D-20.  f/0-t,26.f .60 
(P  & S 4 -*(tsp) 
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Turbulent  C2H4  Jet  Flames 


Fig.  2,  Mean  soot  volume  fraction  aJong  the  axis  compared  to  measurements  by  Neill 
and  Kennedy  (1990)  for  a turbulent  C^//*  flame  with  zero  pressure  gradient. 


348 


Needed  Improvements  of  PDF  Methods  for 
Combustion  Applications 

■ Computationally  Expensive 

■ Direct  Calculation  of  Detailed 
Chemical  Kinetics  - Not  Feasible 

Capabilities  of  Reduced 

Reaction  Mechanisms  ?? 


■ Primitative  Status  of  Mixing  Model 

■ Interactions  between  Turbulence 
and  Chemical  Reactions 
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Development  of  " Rational " Models 
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TURBULENCE  MODELING  FROM  AN  APPLICATIONS 

PERSPECTIVE 
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HOW  SENSITIVE  ARE  RESULTS  TO  BOUNDARY  CONDITION  IMPLEMENTATION 
• e.g.  PRESSURE  b.  C.  ON  A ROTATING  SURFACE 


FOR  TURBULENCE  MODEL  EVALUATION  IT  IS  DIFFICULT  TO  SEPARATE 
NUMERICAL  INACCURACIES  FROM  MODEL  LIMITATIONS 


2 

UJ 

LLi 


o 

h- 

< 

< 

Q 


CD 

CD 

111 

r) 


< 

o 

HI 

z ^ 

£ -j 

QC  LLI 
UJ  Q 


Q 

LU 

> 

o 


LU 
Q 

_J 
LU 
Q 

o 

3uj 


; cc 

s 


CC 

s 

z: 

LU 

2 

LU 

CC 


< 

dO 

Q 

LU 

H 

2 

^j 

1 

2 
Q 
CC 


CL 
O 
oc 

CL 
CL 
< 

I—  X 


</0 


LU  m — 


> 

CD 

LU 

_ z 
X o 

o Q 

LU  £ 
52 
co  < 

. o3 

_ CO  UJ  < 
< UJ  CO  > 
GO  d 13  LU 


UJ 
LU 
Q 
O 

O _ 
LU  O 


LU 

O 

3 

_J 

O 

2 

o 

_J 

3 

O 

X 

to 

2: 

o 

i- 

< 

3 

—I 

< 

> 

LU 


LU 
CD 

UJ 

< 

J- 

z 

UJ 

cc 

UJ 

u. 
u. 

D 

< 

LU 

o 

Q- 


LU 

D 

_J 

< 

> 


cn 

UJ 


UJ 


CO  5 


UJ 


o 

CO  5 
UJ  < 

(-  CO 
n LU 


UJ 


< 

O 

UJ 


5* 

S uT< 
O ^2 
8<3g 

iss 


O 

a. 

o 

-J  CO 

< o 

ii 

^Ul 
UJ  H 

H cc 
2:  < 
O X 
po 

ii 

O £ 
o 00 


LU 

5 

cn 

UJ 

Z W 
ui>^ 

O O p 

cr  < £ 
uj  cc  O 

> 3 UJ 

zom 
o o < 
o < u) 


Z 

LU 

1 

CC 

UJ 

CL 

X 

UJ 

Q 

UJ 

z 

U- 

LU 

a 


LU 


CO 

!£ 

CO 

f 

z 

o 


X 

o 

LL 

o ^ 

lu5 

^3 

z< 

-l> 
< LU 

o_, 

*8 

LUO 

X Cl 
(XI  CL 
X 3 
I-  CO 


9UOOZSU-16-SW 


365 


Rockwell  International 

Rockeldyne  Division 


WHAT  DOES  INDUSTRY  DO? 
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Center  for  Modeling  of  Turbulence  and  Transition 
Workshop  on  Engineering  Turbulence  Modeling  - 1991 


DISCUSSION 


D.  Wilcox  (to  Sindir) 

I understand  that  there  is  some  concerns  with  NASP  contractors  that  RSM 
is  too  expensive.  I find  only  a 20%  increase  in  CPU  time  to  compute  a full 
ESM  with  the  newer  algorithms,  relative  to  the  two  equation  models. 

M.  Sindir  (reply) 

The  model  that  you  came  up  with  is  not  a tool  unless  it  gets  into  the  method- 
ology of  established  codes  that  contractors  use  for  validation.  And  that  is  a 
major  activity. 

T.  Gatski  (to  Sindir) 


You  describe  an  extremely  complex  situation  and  then  use  a Baldwin  Lomax 
model;  I have  a confidence  level  of  zero  in  that  calculation! 


M.  Sindir  (reply) 

I know  and  agree,  but  that  is  what  is  being  done  in  industry.  But  time 
constraints  keep  this  problem  from  being  handled  properly. 

A.  Hsu 

I feel  that  two  items  ought  to  be  added  to  the  list  of  tasks  to  pursue  in 
PDF  modeling.  First,  we  have  to  examine  PDF  for  high  speed  flows,  like 
flows  with  shocks.  Second,  for  the  particle  Monte  Carlo  method,  problems 
involved  in  solution  over  a realistic  geometry  should  be  addressed. 
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Center  for  Modeling  of  Turbulence  and  Transition 
Workshop  on  Engineering  Turbulence  Modeling  - 1991 


A.  Singhal  (to  S.  Pope) 

How  do  we  transition  this  new  and  evolving  technology  into  industry?  In 
the  first  stage  of  transition,  you  introduced  the  joint  PDF  of  velocity  and 
composition,  not  just  the  first  level  of  PDF  which  was  shown  in  the  hierarchy 
to  be  just  composition.  I’m  curious  why. 

S.  Pope  (reply) 

The  reason  for  that  is  that  the  numerical  algorithm  for  joint  PDF  of  velocity 
and  composition  is  really  simple  and  more  economical  than  just  composition 
alone.  It  sounds  strange,  but  the  reason  is  that  in  the  PDF  for  composition, 
the  diffusion  terms,  turbulent  transport, 


Center  for  Modeling  of  Turbulence  and  Transition 
Workshop  on  Engineering  Turbulence  Modeling  - 1991 
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Unconventional  Turbulence  Modeling 
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WORKSHOP  ON  ENGINEERING  TURBULENCE  MODELING  (Aug. /21, 22/1991) 

Center  for  Modeling  of  Turbulence  and  Transition,  ICOMP,  NASA  Lewis 

77^  y 

The  Present  State  of  DIA  Models 

and 

Their  Impact  on  One  Point  Closures 


Akira  Yoshizawa 

Institute  of  Industrial  Science 
University  of  Tokyo 


I.  Objectives 

A,  Outline  of  DIA 

B,  Outline  of  TSDIA  (two-scale  DIA)  and  some 
suggestions  to  turbulence  modeling 

C,  Proposals: 

<D  Helicity  for  the  study  of  the  effects  of 
swi rling  and  cross  flow 
(§  Density  variance  for  the  study  of  highly 
compressed  flows 


PRECEDING  PAGE  BLANK  NOT  FILMED 
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I,  Basic  Laws 
Compressible  fluid 

(3/H)  p + V-  (pu)  = 0 

{d/H)  puj  + (d/hi)  (pujUi)  = - (3/9 Xi)  p 
+ (3/3 xj)  pi s j i 

(3/3 1)  pe  + V"  (pue)  = - pV'u  + $ + V'  (1V0) 
where 

pi  Density;  u:  Velocity;  p:  Pressure; 
e:  Internal  energy;  0:  Temperature; 
p;  Viscosity;  1:  Heat  conductivity; 

Sij  = 9 u j/9 x i I 3 u i/3 x j - (2/3)V,u5iJ 
^ = /» [ (9uj/3xi) 2 + (1/3)  (V-u) 2] 
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MU  iii  mi 


[Note]  Thermodynamic  relations  for  a perfect  gas! 

P = pR0  (R:  Gas  constant) 

= (y  - 1)  e (y:  Ratio  of  specific  heats) 

e = Cvfi  (Cv I Specific  heat  at  constant  volume) 
Incompressible  approximation 


V-u  = 0 

(a/at)  ui  + (a/ax j)  (ujUi)  = - (a/axi)  P + vAui 

where 

p/p  ^ P 

v = p/p:  Kinematic  viscosity 
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II,  Outline  of  DIA 

Premise!  Vanishing  mean  velocity 
Homogeneity 

<iii  (x : t)  Uj  (x'  ;t')>  = <ui  (x  - x' : t)  uj  (0;t')> 
Independence  of  the  coordinate  origin 
^ Infinite  or  periodic  region 
^ Fourier- integral  or  -series  representation 

Fourier  integral 

f (x ; t)  = /f  (k;t)  exp  (-  i k • x)  dk 
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Homogeneous  turbulence 


k*u  (k;  t)  = 0 

(}/it)  Ui  (k;t)  - ikj//8  (k  - p - q)  dpdq 

2 

x iij  (p ; t)  ui  (q ; t)  = ikiP  (k;t)  - vk  Ui 

[Note]  Elimination  of  pressure! 

P (k:  t)  =-  (ki  kj/k2) //S  (k  - p - q)  dpdq 
X ui  (p|t)  uj  (q 1 1) 

[Note]  Green'  s function! 

(a/at)  Gi (k;t,  t')  - iMimn  (k)Jfi  (k  - p 

X um  (pj  t)  Gnj  (q 1 1,  t ) - vk  Gij  (k , t, 
= Dij  (k)  5 (t  - t') 

where 

Du  (k)  = iij  - kikj/k2 

Mi  jk  (k)  = [kjDik  (k)  + ktDi  j (k)  ] /2 


(k|  t) 


q)  dpdq 
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Fundamental  variables 


Qio  (k ; t,  t' ) = <Ui  (k:  t)  Uj  (k'  ; t' ) >/5  (k  + k' ) 

Gu  (k ; t,  t' ) = <0!/  (k; t,  t')> 

Difficulties  in  incorporating  inhomogeneity 

| Necessity  of  the  orthogonal  function  satisfying 
the  noslip  condition 

3 

[Note]  Dannevik'  s work: 

Turbulent  Ray  lei gh-Bena rd  convection  between  two 
parallel  plates  (no  mean  flow) 

(§  Coexistence  of  slow  modes  (mean  field)  and  fast 
modes  (fluctuation) 

— ► Simultaneous  treatment  of  different  modes 
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IV,  Outline  of  TSDIA  (Two ~ S c a 1 e DIA) 

Departute  from  a complete  two-point  scheme 

® Passive!  Difficulty  in  dealing  with  a boundary 

® Positive:  Difficulty  of  obtaining  "formulae" 
applicable  to  general  flows 

Two  scales 

f (=  x),  X (=  5 x)  ; x (=  t),  T (=  St) 
where 

8;  A small  scale-expansion  parameter 
Then 

f = F (X ; T)  + f'  (£,  X;  T,  T)  f F = <f> 
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Two-scale  expressions 


v* *u'  = - 5 a/ih)  u/ , vt  = (a/ati) 

(a/aOui'  + (u-v{)ui'  + (a/afi)p'  - va^' 

= j[-  Uj'  (a/axj)  Ui  - (d/dt)  Ui1  - (a/3Xi)P' 

- (3/aXj)  Uj'  ui'  ] + v — related  terms] 

[Note]  The  effects  of  slow  modes : 

Direct  effects  : Through  U 

Indirect  effects:  Through  X and  T in  u' 

Fourier  representation  of  t 

f'  U.Xjf.T)  = / f'  (k,  X;  I,  T)  exp  [-  ik*  (£  - Ut]  dk 
Scale  expansion 

f'  (k,  X;  T,  T)  = ln=oi  in1  (k,  X;  T,  T) 

[Rote]  Lowest  — order  or  basic  field  uo  * 

The  same  system  of  equations  as  for  homogeneous 
turbulence,  except  the  X and  T dependence 
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Isotropic  and  helical  field 


<uBi'  (k,  X;  t,  T)  uBj'  (k' , X;  t'  , T)  >/S  (k  + k'  J 
= Di  j (k)  Qb  (k,  X;  t,  t' , T) 

+ (i/2)  (km/k2)  s i jmHE  (k,  k' , X;  t,  t' , T) 

_ i 

I — 

helicity  effect 

[Note] 

<uB'  72)  - /Qb  (k,  x;  T,  t,T)dk 

<uB'  ‘Wb1  > = /Hb  (k,  X;  r,  t,  T)  dk,  (Ob'  = V X uBf 

Important  correlation  functions 

Calculation  of  the  Reynolds  stress  etc.  using 
DIA  based  on  Qb  and  Gb 

Extended  eddy- viscosity  representation  for 
the  Reynolds  stress 
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V.  Main  Results  from  TSDIA : No  Helicity 


Reynolds  stress 


- <Ui'  uj'  > = - (2/3)  K« i j + i.Sij 
- 2;  n — 1 f n [Tnij  “ (1/3)  TnmiJ  i j] 


U (D/D  t)  Si  j + 


where 

K = <u'  Z/2> 

Sij  = dUj/dxi  + dUi /hi 

Tiij  = (iDi/axJ  (dUj/dxJ 

T2ij  = [(dUi/dxJ  (dlL/dxj) 

I (dUj/dxm)  (dUm/dxi)]/2 
Tsij  = (3Um/dxi)  (dUm/dxj) 

5 

[Note]  See  Speziale  and  Rubinstein  and  Barton 


i 

R 

i 
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Turbulent  scalar  flux 

<ux'  0'  > = - ti  j (i/hi)  0 
where 

ku  = icS i j - [y i (3Ud/3xi  + dUi/dxj) 

- y 2 (3Uj/xi  - dUi/xj)  ] 

8 

[Note]  See  Rubinstein  and  Barton 
Some  other  suggestions 

. __  9 

($  Triple  correlations 

<(u'  2/2)  u'  > = - vkVK  + »K'Vt 

® Equations  for  the  dissipation  rates  of 

7 

energy  and  scalar  variance  {i  and  t e) 
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W.  Proposal;  Helicity 
Motivation 


Explanation  of  the  generation  of  'white  spot'  in 
Saturn  (the  spiral  vortical  structures  along 
the  equator) 

Fig.  1;  Saturn's  huge  white  spot 

(NASA  Hubble  Space  Telescope) 


What  is  the  helicity? 

(3/3 1)  u = u X « - V (p  + u /2)  + vAu 

u|{i)-^uX(1)  = 0 (no  energy  cascade) 

Helicity  u*(o;  A measure  of  the  break  of 

reflectional  symmetry  in  flow 


Fig.  2.  Helicity 
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Importance  of  helicity 


% A measure  of  'smallness'  of  energy  cascade 

§ A 'conserved"  quantity  in  the  absence  of  the 
injection  and  loss  due  to  viscous  effects 

[Note]  <u*(i)>  = U*Q  + <u'  * > 

% No  mean  helicity!  U • Q - 0 (U  1 fi) 
Two-dimensional  mean  flow 
channel  flow,  jet,  wake,  mixing  layer,  etc. 

§ Finite  mean  helicity!  U*Q  + 0 

Swirling  flow,  Three-dimendsional  mean  flow 

[Note]  Swirling  flow! 

Question!  Why  do  the  eddy-viscosity  models 

break  the  swirling  motion  so  fast? 

Answer  ! No  consideration  of  helicity  or  the 

decreasing  effect  of  cascade  (virtual 
decrease  of  eddy  viscosity) 
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[Note]  Cross-flow  effects! 

§ Lag  of  the  turbulent  stresse  in  response  to 

the  cross-flow  gradient  (also  important  in 

11,  12 

aerodynamical  flows  ’ ) 


(§  Curvature  effects  leading  to  secondary  flows 

13 

Three-equation  model 


K = <u'  /2>:  Turbulent  kinetic  energy 
i : Energy  dissipation  rate 

H = <u'  •(/>:  Turbulent  helicity 
[Mean  equation] 

V-U  = 0 


(D/D t)  Ui  = - (d/d x i)  P i (d/dxjjRji 


[Reynolds  stress] 


Rji  = - Cui'  u/> 

= - (2/3)  K« ij  + ve  (3Uj/ixi  + JUi/Jxj) 
- [Qi r j + Q j y i - (2/3)  Q • r 5 ± j] 

where 

Ve  = Cxi  2A,  r = c2(k4/e3)VH 


[Note]  Symmetry-breaking  factor: 

Reynolds  stress:  Ref  1 ect i ona 1 ly  symmetric 

Necessity  of  another  symmetry-breaking 
factor,  that  is,  inhomogeneity  such  as  VH 
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[K  equation] 

(D/D t)  K = PK  - £ + V*  [C3  (K2/0  VK] 
where 

Pk  = Rij  (dUj/Jxi) 

[H  equation] 

(D/D t)  H = Ph  - £h  + V-Th 
where 

Ph  = Rij  (8Qj/3xi)  - fii  (3/xj)  Rji 

£ h = C4(£/K)H 

Th  = Kfi  + Cs(K2/£)VH 

[f  equation] 

(D/Dt)  £ = C6  (c/K)  Pk  - C7  (£2/K)  + C8  (K2/£)  Q‘VH 
+ V’  [Cg  (R2/ £ ) V £ ] 
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Swirling  flow  in  a pipe  [ (r,  6,  z)  5 z‘.  axial] 

- <ur'  u/>  = ver  (3/a r)  (U,/r)  - lQ.iH/3r 

- <(u r'  Hz'  y - v e (3Uz/3r)  - lUz^H/^r 

- <u g'  .:>  = Ve  (3Ue/3z)  - ?!UH/3z 
where 

4,3.. 

i}  = CiK  A > 0 

1 4 

Comparison  with  observation 

Fig.  4.  Mean  velocity 
(Broken  lines:  weak  swirl 

- <u r'  us'  > = ver  (3/3 r)  (Uo/r)  - 1 Q o 3 H/ 9 r 

A B 

Around  r - 0.  7 

A > 0 

B>0  (Qo  = - dUz/d  r < 0,  d H/  9 r < 0 ) 

[Note]  H = 0 at  the  wall 
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Fig.  5.  - (u/  Ufl'  y 


region) 


3H/9  r < 0 near  it 


Application  to  SGS  modeling 

Appearance  of  steamwise  vortices  (streaks) 
dependent  on  the  strength  of  shear 

— ► Nonvanishing  helicity  on  the  SGS,  hut  not  ; 

on  the  ensemble  or  time  mean  ; 

•+■  Importance  in  the  SGS  modeling  \ 

Incorporation  of  the  helicity  effect 

-4  Virtual  change  of  the  Smagorinsky  constant 
in  channel  flow,  mixing  layer,  isotropic 
flow  in  accordance  with  the  strength  of  j 

shear 
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in  i mi  mi  nun  iui  mi  him  mi  mi  iii  iiiiiiiiiii  mil  mmwwi 


W,  Proposal'.  Density  Variance 

Prominent  difference  between  incompressible  and 
compressible  turbulence 


Decelerated  streamwise  velocity  effects 

(§)  Compressible  (shock  wave)  . 

Decrease  in  turbulence  level  ; 

Increase  in  temperature  etc,  : 

§ Incompressible:  j 

Increase  in  turbulence  level — 

Mass-weighted  mean 

f = {f}  = <pf>/p“»  7 = 

{'  = f - f 
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[Mean  equation] 


(a/at)  ? + v-  (pS)  = o 

(a/at)pui  + (a/axj)  pu'jUi  = - (a/xi)  [(r  - i)pe] 
+ (a/du)Rji 

(a/at)  pe  + V*  (pue)  = - (r  - 1)  peV*u  + p"(f  + l) 

+ V-  (-  H) 

where 

Rs  j = - <puj'  uj'  > = - J {u 3 ' uj' ) 

H = <p  u'  e'  > = p ju'  e' } 

z = - (r  - lXpe'V-u 'Wp  = - (y  - D (e'  V-u' } 

RijI  Reynolds  stress 

H : Internal  energy  flux 

l I Fluctuating  dilatation  effect 
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[Note]  Importance  of  V 

Trace  of  the  pressure-strain  correlation 
(no  contribution  in  the  incompressible  case) 

[Note]  Inference  of  l\ 

Large  positive  in  a highly  compressed  region 
Virtual  increase  of  energy  dissipation 

j i 16’  17 

Three-equation  model 

K = {u'2/2},  f,  Kd  = <p'  2>:  Density  variance 
[Reynolds  stress] 

Rij  = - (Z/8)  pKiij  + P V e S i j j Ve  = CiK2/f 
[Internal  energy  flux] 

H = - p/CeVe,  Ke  - C2K  /t 
[Fluctuating  dilatation  effect] 

l = C3  (f/K)  (Kd/p2)  e 
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[K  equation] 


(a/a t) k + v-  (puK)  = pk  - p (t  + z)  + v-Tk 

[f  equation] 

(a/at)  t + V-  Rut)  = C4  (t/K)  Pk  - CsHt/K)  (t  + z) 
+ V'T, 

[Kd  equation] 

(d/at)  Kd  I V-  (uKd)  = - KdV’u  - CG(f/K)Kd  + V-Td 

i i 

exact ! 

Can  the  turbulence  level  decrease  behind  a shock 
wave? 

V*u  < 0 near  a shock  wave 
"►  Production  of  Kd  (Kd  equation) 

"►  Larger  e,  larger  l (than  elsewhere) 

Decrease  in  K and  1 (K  and  t equations) 
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428 
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Time  evolution  o(  the  turbulent  kinetic  energy  In  homogeneous  shear  dovr 
Relaxation  model;  o Large-eddy  simulation  of  Bardina  et al  for  eo/SK^O.296 

Figure  2 
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BACKWARD-FACING  STEP: 


3 


Y/H 


X/H 


(a)  Streamlines 


(b)  Dimensionless  mean  velocity  profile 

- Compulations  with  Isotropic  eddy  viscosity; 
o Experiments  of  Kim  et  zl,  1980;  Eaton  & Johnston,  1981) 

Computed  mean  flowlleld  for  the  new  RNG  K-z  model 
[E  = 1:3;  Re  = 132,000;  200x100  mesh] 

Figure  4 
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BACKWARD-FACING  STEP: 


(a)  Turbulence  Intensity 


(b)  Turbulence  shear  stress 

ComputedturbulencestressesforthenewRNG/C-emodel[E=1:3;Re  = 132,000; 

200x100  mesh;  computations  with  Isotropic  eddy  viscosity; 

o experiments  of  Kim  ef  al,  1980;  Eaton  & Johnston,  1981] 

Figure  5 
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BACKWARD-FACING  STEP: 


X/H 

(0)  Streamlines 


(b)  Dimensionless  mean  velocity  profile 

( Computations  with  anisotropic  eddy  viscosity; 

o Experiments  of  Kim  et  at,  1980;  Eaton  & Johnston,  1981) 


Computed  mean  flowfleld  for  the  new  RNG  K-e.  model 
[E  = 1:3;  Re  a 132,000;  200x100  mesh] 

Figure  6 
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BACKWARD-FACING  STEP: 


X/H  : 1,33  5.33  7.67  8.55  10.33  15.66 


(a)  Turbulence  Intensity 


0.0  0.010.02  0.010.02  0.010.02  0.01  0.02  0.01  0.02  0.010.02 


(b)  Turbulence  shear  stress 

ComputedlurbulencestressesforthenewnNG/C-r.model[E=:1:3;Re  = 132,000; 

200x100  mesh;  compulations  with  anisotropic  eddy  viscosity; 

o experiments  of  Kim  ef  af,  1980;  Eaton  & Johnston,  1981] 

Figure  7 

433 


- ft  £-*1 

N92-24535 


H 

Q 

CC 

O 


if) 


u. 

O 

LU 
h- 
< 
H 
C 0 

H 

Z 

LU 

if) 

LU 

CC 

Q. 

LU 


if) 

£ 

O 

IL 

>■ 

CC 

Q 

LU 

< 

H 

LU 

Z 

h 

LU 

if) 

7 

5^0 

3 

U. 

< i 5= 
CD  O S 

* ic  O 

O 

z < q' 

0 

z 

< LU  7 

2 </>  5 

. LU  5 

_J 

CC  DC  UJ 

LU 

< CD  > 

O 

O 

2 

o > y 

guJO 

UJ 

< 

o 

if) 

Z 

< 

LU 

Z 

J 

3 

CD 

CC 

3 

1- 

435 


PRECEDING  PA3E  BLANK  NOT  FILMED 
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THE  DISSIPATION  TRANSPORT  EQUATION  SHOULD  NOT 
BE  BASED  ON  LOCAL-EQUILIBRIUM  HYPOTHESIS  BUT  ON 
RAPID-DISTORTION  HYPOTHESIS.  NO  SUCH 
DIFFERENTIAL  EQUATION  IS  CURRENTLY  AVAILABLE. 
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REPRESENTS  THE  INFLUENCE  OF  THE  MEAN-STRAIN- 
RATE  DEVELOPMENT  ON  AND  THE  THIRD 

ACCOUNTS  FOR  ANY  IMPOSED  LENGTH  SCALE. 
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(Tradu  et  al . 1991) 


Effect  of  Frequency  on  Turbulence  Intensity  - From 

Mizushina  et  al.  (1973) 
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CD-90-50119 


Aiodu/ac/on  of  wall  shear  stress 


(Ronneberger  & Binder  H91) 
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The  Calculated  Cycle-Variations  in  the  Reynolds  Shear 
Stress  at  Re  = 50,000  in  Comparison  with  Tu  and 
Ramaprian's  (1983)  Data  for  a Pipe  Flow  at  y/R  = 0.035 
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TO  OBTAIN  DEFINITE  PREDICTIONS  IT  IS 
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EQUATION. 


Piosent  rapid-distortion  model 
Quasi-steady  K-t.  model 
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Man  U llanratUy  (1906),  Q<  ar»'  Ramaprian  & Tu  (1903),  D-  (a)  Amplitude 
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THE  ROLE  OF  EXPERIMENTS 
IN 

TURBULENCE  MODELLING 


by 


William  K,  George 
Turbulence  Research  Laboratory 
State  University  of  New  York  at  Buffalo 


PRECEDING  PAGE  blank  NOT  FILMED463 


1. 

2. 

& 

Some 


THE 


recent  references 


w.  K.  George 


W.  K.  George 
D . B . Taulbee 


"Governing  Equations,  Experiments, 
Experimental  Fluid  Thermal  Science 


and  the  Experimentalist 

- 1990 


"Designing  Experiments  to  Test  Closure  Hypothesis" 
to  ap^ar8Exp?mental  Fluid  Thermal  Science  - W* 
(Also  Proceedings  Synmposium  on  Tur  u 
Dubrovnik  1990) . 


fvppri mental  PiletTmjas  --  just  a few  of  many 
. the  Round  Jet  - what  is  its  growth  rate? 
the  Turbulent  Boundary  Layer 
the  Plane  Wake 
. the  Dissipation 


ROUND  JET 


6-t/o  “ 0.095  X 
not  0.085  X 

Illustrates  one  of  the  biggest  problems 


for  modeller 


Who  to  believe? 


Problem  for  experimentalists  is 
- How  to  overturn  old  results? 

- unpublished  since  disagree  with  earlier  results 

- unaccepted  since  unpublished 


— unused  since  unknown 


. Problem  for  everyone  - Isn't  there  better  way 

— to  get  results  out 

— to  purge  old  results 


j 


m 

! 


464 


LOCAL  ISOTROPY  ??? 


- Local  Axisymmetry? 

- Equl -partition  of  Dissipation? 


465 


0314  P06 


L.  W.  B.  Brown*.  Jt.  A.  Antonia  and  D.  A . Shah 


».< 


I.:’ 


e^KlO-* 


r 


mi 


^0 


Piorti  4.  Distribution  acror*  the  wake  of(fl/f/dy)* . the  r.m.i.  of  derivative*  of  velocity  fluctuation*. 
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O.ft  — 1 (figure  Id);  □,/  ■ to  (figure  1c);  ^ « a.  u»ing  two  single  wire*. 
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THE  PIANE  WAKE 

- there  is  no  wake  which  is  independent  of  generator 


iiiiiiiiiilllinll  i ii  i !M 


M 


Ftfer*  5. 


Variation  cootorllM  velocity  deficit  and  ^ilf-widtfc 
w1tJ»  distinct  for  t hrtt  wak*  generators;  a .airfoil; 
a,  70X  solidity  screen;  0 . wild  strip  (fro* 
Wygnanskl  ft  al.  1SS6). 


Figere  6. 


ynolds  stress  nomalized  by  centerlliM  velocity 
r Dm  solid  strip  and  airfoil  (Vygnaaskl  et  al.  1906). 
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Plgvrt  7.  Normalized  tirfevliact  latmHjf  profiles  for  the  three 
generators  of  Figaro  S Qfrgainshl  ot  el.  1991). 

Tns  Axzsymmctiiic  Uakc:  A Flow  Wkicw  Docs  Not  Evolvc  at 

Constant  Reynolds  Nunsck 


Tbo  axlsymetrlc  wake  pro  loots  on  Interesting  contrast  to  the 
axlsymetrlc  jot  and  piano  wake  flows  described  above  In  that  It  does 
not  evolve  at  constant  Reynolds  number  (as  will  be  seen).  As  x 
consequence,  the  nature  of  the  assumptions  regarding  the  dissipation 
will  be  seen  to  predict  two  quite  different  asynptotlc  developments. 

There  appears  to  be  experimental  evidence  for  both  forms  In  different 
experiments,  which  raises  an  Interesting  question  as  to  how  the  flow 
chooses  one  form  or  another.  An  Interesting  possibility  is  that  the  1 

flow  evolves  from  one  state  to  another  as  the  Reynolds  nuaber  changes. 

These  possibilities  will  be  discussed  In  more  detail  below  following  the 
analysis. 

The  equations  of  motion  describing  the  axlsymetrlc  wake  to  first 
order  can  be  shown  to  reduce  to. 


u.  If  («-«.]-- £$?  «■  « (78) 

where  U.  Is  the  undisturbed  speed  of  the  free  stream.  This  can  be 
integrated  across  the. flow,  iojrltld  the  Integral  constraint. 


**  (IHJJrdr  - (79) 

/o 

where  # Is  defined  to  be  the  omentum  tkickaess. 
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This  persistent  effect  of  initial  conditions  may  be  widespread! 


471 


The  'Up! 1 -established*  Turbulent  Boundary  Layer 
. Is  there  really  a log  layer? 

Should  comparisons  be  made  with  these  correlations? 
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Fig.  1.1.1  On  the  determination  of  the  friction  coefficient. 
from  KLINE  ct  al.  (1967). 


Mjx.  rf  -fo  Ca/c*,/^  yci/  re. 

/A  ^ — *** 
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Fig.  1.2  Velocity  derivative  from  Direct  Simulation: 
from  SPALART  (1988) 
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Fig.  4.4.3. 1.1  Shape  factor. 

Smith  and  Walker,  x = 15.75  in,  27.75  in,  39.75  in.  51.75  in. 


480 


.iiiNiiinm 


Fig.  4.4.6. 1.1  Friction  coefficient  as  a function  of  x. 
Smith  and  Walker  (53  points). 
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The  Past 


Without  experiments,  there  would  be  no  single -point  turbulence 
models . 


DNS,  LES,  etc.  have  helped  augment  data  base  in  recent  years. 
Contributions  have  been  important,  but  relatively  minor  when  taken 
as  part  of  the  whole* 


The  Future 

. Question:  Is  this  balance  likely  to  change? 

. No  question  simulations  will  play  increasingly  important  role 

— ability  to  produce  difficult  quantities  (for  expts) 

(pVj  , pressure  strainrate,  dissipation,  etc.) 

— increasing  Reynolds  number  (still  small) 

— National  economic  agenda  of  U.S.  - we  will  do  what  is 
necessary  to  maintain  market  position 

— indirect  subsidy  of  CFD 

. What  role  will  this  leave  for  expts? 

Experiments  of  value  to  modelling  community  are  very  difficult  to  do 
Most  efforts,  however  well- intended,  do  not  measure  up! 

Reasons  complicated  but. . . high  on  the  list. . . . 

, Inexperience  of  investigator. 

. Ignorance  of  goals. 

Money  usually  runs  out  before  experimenter  learns  to  do  it 
right. 

. These  will  be  Increasing  problems  in  next  decade. 

Biggest  Problem  for  Experimentalist: 

Few  sponsors  have  patience  to  see  an  experiment  through  to  its 
completion  - none  in  my  experience! 

As  complexity  of  flows  to  be  investigated  increases  (Bushnell 
comment),  this  problem  will  be  exacerbated! 
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FINAL  THOUGHT 

However  successful  our  closure  efforts  may  be.  we  are  not  solving  the 
turbulence  problem  - we  are  only  being  responsible  engineers. 


We  must  therefore  be  careful  not  to  demean  the  efforts  of  those  who  are 
trying  to  grapple  with  the  real  turbulence  problem  “there  may  be  an 
^erfying  physical  principle  which  will  very  much  affect  what  we  do! 
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NUMERICAL  SIMULATION: 

ITS  CONTRIBUTIONS  TO  TURBULENCE  MODELING 
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Homogeneous  turbulence 
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MODEL  EQUATIONS 
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4-  C^Vt  {dyU) 
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llUil  Kill.  IWiHKittllilil  , ii  ill 


variances 

k and  v2 


493 


Fig.  lb.  k (solid  line)  and  v3  (dashed  line)  from  the  model 
compared  to  DNS  data  (symbols). 


COMPRESSIBLE  TURBULENCE 
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Boundary  layer  questionable 


Fig.  1 Pressure  dilatation  durin 


STOCHASTIC  SHOCKLET  DISSIPATION  MODEL. 


FIGURE  1,  (a)  Sketch  of  shock-like  structure  in  a turbulent  eddy;  (b)  normal 

shock  relations. 


Average  dilatation  dissipation: 


ed<Xi[wL  i^1)  P(m,K)dm] 


and 


€d  — esCdF(Mt,K) 


Total  dissipation  (for  given  kurtosis  K of  m) 


ttot  = es{  1 + cdFk(Mt)} 
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Centerline  Turbulence  Intensities  in  Mixing  Layer 


energy  to  the  normal  shock 
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Stratified  flows 


DIRECT  SIMULATION 
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FLOWS  SIMULATED-CONTINUED 
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Inhomogeneous  Incompressible  Turbulence 
Wall-bounded  flows 

Channel  Curved  channel  Boundaiy  Layer 
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EXAMPLES-MODEL  TESTING/VALIDATION 
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EXAMPLES-MODEL  DEVELOPMENT 


<D 

O 

C 

<D 


cc 


a> 

o 

c 

0) 

D 

JQ 

i_ 

3 

0 

n 

w 

w 

(0 

CD 

<1) 

>» 

CQ 

Q. 

E 

o 

O 

I 

606 


Flamelet  models 


ROTATING  TURBULENCE 

Isotropic  turbulence  subject  to  rotation 
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RESULTS-ROTATING  TURBULENCE 
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ANISOTROPY  MEASURE 
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WALL  BOUNDED  FLOWS 

Near-wall  turbulence  difficult 


Bll 


Good  predictions  of  boundary  layers 


Have  modelled  blocking  effect  with  new  k — e — v model. 
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Center  for  Modeling  of  Turbulence  and  Jra^3ltton 
Workshop  on  Engineering  Turbulence  Modeling  - 1991 


DISCUSSION 


B.E.  Launder  (to  S.  Orszag) 

You  showed  us  values  of  C(i  and  U(2  which  would  give  too  high  a decay  rate 
of  grid  turbulence. 

S.  Orszag  (reply) 

There  is  some  question  about  what  the  decay  rate  really  is. 

B.E.  Launder  (to  S.  Orszag) 

You  talked  about  C<2  but  you  didn’t  say  how  C{1  emerged. 

S.  Orszag  (reply) 

That’s  the  same  calculation. 

S.  Pope  (to  S.  Orszag) 

Calculations  you  did  with  backward  facing  step,  what  boundary  conditions 
did  you  use? 

S.  Orzsag  (reply) 

That  was  not  a full  RNG  calculation.  It  should  have  been  done  using  in- 
terpolation formulation  for  the  various  constants  all  the  way  to  the  wall. 
Instead  it  was  done  using  the  a fit  right  to  the  log  layer. 
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Center  for  Modeling  of  Turbulence  and  Transition 
Workshop  on  Engineering  Turbulence  Modeling  - 1991 


T.  Gatski  (to  W.K.  George) 

Two  things  modelers  are  looking  for  validation  and  calibration.  Bill  refered 
to  the  kind  of  experiments  we  use  for  validation.  The  kind  of  work  you  (J.H. 
Ferziger)  do  with  DNS  has  been  building  block  for  calibration.  Do  you  want 
to  design  the  experiments  to  validate  our  model  or  calibrate? 


W.K.  George  (reply) 


I would  like  to  design  experiments  which  would  invalidate  your  models. 


T.  Gatski  (to  W.K.  George) 


But  that  would  be  destructive  for  both  of  us. 
P.  Spalart  (to  J.H.  Ferziger) 


I differ  with  your  description  of  DNS  as  exact  solution.  I would  like  to  say 
that  my  solution  are  not  exact.  I spend  quite  a lot  of  time  thinking  how  I 
can  keep  the  error  small.  If  I can  double  the  number  of  the  grid  points  in 
each  direction  I would  sleep  better  but  it  would  take  fifty  years  instead  of 
two  years  to  finish  the  simulation. 

J.H.  Ferziger  (reply) 

I should  have  said  that  in  my  talk  that  any  numerical  calculation  is  approx- 
imate and  I hope  we  work  hard  to  keep  them  small.  There  are  errors  due  to 
numerical  methods  and  errors  due  to  the  fact  that  we  have  limited  computer 
time. 
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Center  for  Modeling  of  Turbulence  andTrja?siU°"(,a, 
Workshop  on  Engineering  Turbulence  Modeling  - 19»1 


p.  Spalart  (to  W.K.  George) 


A comment  on  BiU’s  theory  that  it’s  a power  law  instead  of  log  law.  I think 
in  your  original  APS  abstract  you  make  it  sound  like  it  s just  a matter  of 
taste  if  you  use  defect  law  or  your  theory. 


W.K.  George  (reply) 


In  the  original  APS  abstract  I wasn’t  clear  why  the  theory  comes  apart 
asymptotically  and  I feel  much  confident  now  that  the  existing  theory  is 

wrong. 


P.  Spalart  (to  W.K.  George) 


There  is  very  different  Galelian  invariance  to  those  two  theories  and  its  not 
a matter  of  test.  Defect  law  says  we  are  coming  from  the  free  stream  and 
we  don’t  know  how  fast  the  wall  is  moving;  it  may  be  a moving  belt.  lour 
theory  doesn’t  do  that. 


W.K.  George  (reply) 

That’s  right.  There  is  no  question  that  there  is  a lot  of  sorting  out  to  be 
done. 


B.E.  Launder  (to  W.K.  George) 

It’s  very  interesting  that  he  (or  it)  brings  out  into  question  the  universality. 
If  you  got  flows  to  decrease  rapidly  with  distance  from  the  wall  as  you  do 
in  low  Re  channel  flows  there  is  data  going  back  to  fifties  that  your  og-  og 
constant  goes  up.  So  logically  you  would  expect  log-log  constant  would  go 
down  in  adverse  pressure  gradients.  The  implicit  faith  shown  in  sectors  of 
fluid  community  in  the  universality  of  the  log  law  I think  is  misplace  . 


515 


Center  for  Modeling  of  Turbulence  and  Transition 
Workshop  on  Engineering  Turbulence  Modeling  - 1991 


There  is  now  emerged  which  I think  an  excellent  paper  by  Nagano  (at  the 
upcoming  SFC  in  Munich)  showing  what  seems  to  me  a cleax  dependence  of 
log-log  constant  on  shear  stress  gradient.  In  adverse  pressure  gradients  lower 
log-log  constants  than  you  find  in  zero  pressure  gradients.  Bill  gets  unhappy 
at  unacknowledged  at  his  private  discoveries  as  all  of  us  do,  I suggested  this 
in  a paper  about  eight  years  ago. 


W.K.  George  (reply) 


I presented  this  in  1978. 

P.  Spalart  (to  B.E.  Launder) 

I have  results  that  show  that  in  moderate  pressure  gradients  log-log  going 
down.  At  = 50  it  goes  down  by  almost  one  wall  unit  at  /?  = 2 which  is 
not  very  strong  at  all. 

S.  Pope  (to  P.  Spalart) 

You  very  quickly  mentioned  that  you  use  DNS  data  for  guidance  and  not 
calibration.  Could  you  expand  on  that? 

P.  Spalart  (reply) 

If  I calibrate  turbulence  model  for  boundary  layer  based  on  flat  plate  results 
I’ll  get  too  high  Re,  so  I don’t  have  DNS  results  which  I’ll  trust  within  10% 
to  extrapolate  to  high  Re. 

G.  Hwang  (to  J.H.  Ferziger) 

You  talked  about  P.  Durbin  model  that  uses  v2  as  damping  function.  But 
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Center  for  Modeling  of  Turbulence 

Workshop  on  Engineering  Turbulence  Modeling  - 


problem  is  we 
talking  about 
stress  models. 


want  to  use  this  model  in  multi-dimensions.  Then  you  are 
k -v2  and  e equation  which  is  as  complicated  as  Reynolds 
I am  not  disputing  the  model  there  axe  some  good  points  to 


it. 


A.  Yoshizawa  (to  J.H.  Ferziger) 

You  pointed  out  the  possibility  of  LES  as  engineering  tool  but  our  experience 
shows  we  can  not  perform  LES  with  Smorginsky  constant  fixed.  LES  criti- 
cally depends  on  Smorginsky  constant  e.g.  we  perform  LES  o anne 
usine  Ca  = 3.1  but  using  this  constant  we  can  not  simulate  e.g.  backward 
facing  step.  My  opinion  is  without  overcoming  this  difficulty  LES  can  no 
become  engineering  tool. 


J.H.  Ferziger  (reply) 


I aeree  with  you.  There  is  a new  model  which  does  overcome  some  of  these 
difficulties.  I don’t  want  to  say  that  all  the  difficulties  are  overcome  but  we 

have  hope. 


J.  Bardina  (to  A.  Yoshizawa) 

We  have  investigated  DIA,  RDT  and  LES.  LES  is  much  simpler  than  DIA 
You  have  difficulty  with  Smorginsky  constant  but  we  know  that  this  constan 
is  not  right.  It’s  not  universal  for  all  flows  because  only  thing  it  s doing  more 
is  dissipating  more  energy;  that’s  all.  At  high  Re  you  are  putting  more 
energy  at  small  scales  and  there  are  many  other  effects  which  you  have  to 

put  there. 
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Center  for  Modeling  of  Turbulence  and  Transition 
Workshop  on  Engineering  Turbulence  Modeling  - 1991 


E.  Reshot  ko 

Bill  brought  up  many  of  the  problems  in  looking  at  experiments.  There  are 
many  more  which  try  to  fulfill  some  of  the  desires  and  expectations  that  have 
been  brought  up  here  at  the  workshop.  First  any  self-respecting 

experimentalist  is  not  there  to  design  experiments  to  validate  or  calibrate  a 
theory.  First  of  all  he  is  there  to  discover  new  physics.  May  be  a theoretician 
wants  to  see  if  the  physics  is  reproduced  by  the  model. 

But  when  it  comes  to  doing  that  particularly  in  measuring  turbulence  we 
come  up  with  the  problem  of  how  to  measure  at  a point.  We  have  all  these 
wonderful  things  at  a point  when  we  have  probes  which  are  not  a point. 
Recently  we  had  experience  with  multiple  wire  probe  that  showed  our  probe 
was  not  measuring  at  a point  although  our  probe  was  less  than  0.1mm  in 
overall  size.  And  this  problem  becomes  worse  if  one  goes  to  high  speeds. 
I understood  just  a few  years  ago  why  all  good  turbulence  measurers  were 
working  in  large  facilities  and  in  low  speeds  because  only  in  that  way  you  can 
feel  reasonably  secure  that  in  terms  of  wall  units  you  are  operating  at  a point. 
We  tried  running  some  experiments  at  100  ft/sec  and  found  that  our  probe 
was  100  wall  units  which  typically  a spanwise  streak  size.  In  compressible 
flows,  aside  from  increased  speeds  and  increased  probe  dimensions  in  wall 
units,  we  also  have  the  problem  of  calibration  in  transonic  regimes.  It’s 
not  that  we  don’t  know  how  a hot-wire  works  in  transonic  regimes.  It’s  so 
sensitive  to  Mach  number  in  transonic  regime  that  there  is  not  a way  of 
saying  it’s  reliable.  I am  worried  about  the  double  and  triple  correlations 
in  boundary  layers  with  the  present  probes.  One  of  the  things  H.  Nagib 
is  doing  is  looking  at  probe  miniaturization  and  I encourage  this  but  until 
then  I think  prospects  of  getting  detailed  compressible  flow  measurements 
are  dim. 
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J.H.  Ferziger  (to  all  participants) 

I thought  it  would  be  interesting  to  throw  out  at  the  modelers  that  what  do 
they  think  is  missing  in  the  experiments,  simulations  and  theory? 


B.E.  Launder 

I ponder  from  time  to  time  about  these  nice  homogeneous  flows  which  people 
use  to  come  up  with  constants  in  dissipation  equation  as  Steve  Orzsag  was 
talking  about.  Question  comes  to  my  mind  that  its  the  variation  in  mho- 
mogeniety  which  we  are  interested  in  looking  at.  The  variation  m sPatia 
length  scales  ought  to  enter  in  our  closures  in  ways  other  than  the  diffusive 
like  terms.  That  is  to  say  perhaps  if  we  are  thinking  of  dissipation  equation 
having  adjacent  to  one  another  layers  of  different  length  scales  are  going  to 
be  promoting  spectral  transport  of  energy  removed  more  readily  than  you 
find  in  homogeneous  flows.  So  I ask  myself  if  DNS  can  help  clarify  this. 


J.H.  Ferziger  (to  B.E.  Launder) 

Are  those  relatively  simple  inhomogeneous  flows  in  that  regard. 


B.E.  Launder  (reply) 

I think  simpler  inhomogeneous  flow  you  are  talking  about  is  channel  flow 
where  everything  is  so  dominated  by  the  fascinating  structure  of  the  near 
wall  sublayer.  If  you  could  do  simulation  away  from  wall  where  low  Re 
dissipation  issues  are  dominated. 


J.H.  Ferziger  (to  B.E.  Launder) 


Then  there  are  mixing  layers. 
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B.E.  Launder  (reply) 

Maybe  mixing  layer  results  will  be  valuable.  I don’t  think  they  ought  to  be- 
long to  this  question.  That  is  your  dissipation  equation  needs  to  be  different 
in  inhomogeneous  from  that  in  homogeneous. 

J.H.  Ferziger  (to  B.E.  Launder) 

We  tried  doing  flow  simulation  of  experiments  of  Warhaft.  May  be  we  can 
collaborate. 

B.E.  Launder  (reply) 

Maybe  those  experiments  themselves  will  answer.  There  are  no  mean  velocity 
gradients  in  that  experiment  John  (to  J.L.  Lumley). 

J.L.  Lumley  (to  B.E.  Launder) 

No  mean  velocity  gradients  - just  a gradient  of  scale.  I was  just  going  to 
draw  attention  to  that. 

J.  Bardina  (to  B.E.  Launder) 

Brian  are  you  suggesting  that  homogeneous  flows  axe  not  a valid  test  for  k — e 
model  since  it  would  not  account  for  inhomogeneous  part  in  shear  flows. 
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B.E.  Launder  (reply) 

We  should  look  at  homogeneous  flows  as  building  blocks  and  may  bewe 
shouldn’t  say  going  from  homogeneous  to  inhomogeneous  flows  just  adds  a 
diffusion  like  transport  term  but  may  be  adds  other  as  well. 


J.  Bardina  (to  B.E.  Launder) 


The  k - e model  worked  good  for  homogeneous  flows  we  tested  and  it  wasn  t 
tuned  for  these.  We  looked  at  homogeneous  shear,  and  plane  strain  flows 
and  it  did  very  well.  It  didn’t  do  well  for  rotational  flows  because  effect  of 
rotation  isn’t  accounted  for. 


W.K.  George  (to  B.E.  Launder) 

I think  you  are  right  Brian.  If  you  look  at  Antonio’s  dissipation  results  all 
but  one  derivative  is  way  out  of  line.  I have  come  to  believe  that  it  associated 
with  inhomogeniety.  And  if  you  look  at  each  term  in  the  equation,  1 t e o 
is  truly  locally  homogeneous  you  can  not  produce  any  of  those.  This  question 
is  best  resolved  by  DNS  of  inhomogeneous  situation  at  low  Re. 


J.  Weinstock  (to  B.E.  Launder) 

I wouldn’t  have  any  doubt  that  inhomogeniety  would  cause  changes  other 
then  diffusive  transport.  Nature  of  the  change  is  such  that  we  can  not  tell 
until  we  do  it.  I did  a calculation  where  I accounted  for  strong  time  variation 
where  the  turbulence  energy  is  changing  at  times  of  order  of  eddy  circulation 
time  and  a simple  result  came  out  of  that.  It  wasn’t  diffusive  transport^  Eddy 
circulation  becomes  function  of  turbulence  time  scale.  So  the  coefficients 
involving  damping  become  function  of  rate  of  change  of  turbulence.  t in 
if  we  put  spatial  inhomogeniety  we  would  come  up  with  something  related 
to  that.  It’s  obviously  not  diffusive. 
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